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Abstract.  Biomass  burning  (BB)  is  a  major  global  source 
of  trace  gases  and  particles.  Accurately  representing  the  pro¬ 
duction  and  evolution  of  these  emissions  is  an  important  goal 
for  atmospheric  chemical  transport  models.  We  measured  a 
suite  of  gases  and  aerosols  emitted  from  an  8 1  hectare  pre¬ 
scribed  fire  in  chaparral  fuels  on  the  central  coast  of  Cali¬ 
fornia,  US  on  17  November  2009.  We  also  measured  physi¬ 
cal  and  chemical  changes  that  occurred  in  the  isolated  down¬ 
wind  plume  in  the  first  ~4  h  after  emission.  The  measure¬ 
ments  were  carried  out  onboard  a  Twin  Otter  aircraft  outfit¬ 
ted  with  an  airborne  Fourier  transform  infrared  spectrome¬ 
ter  (AFTIR),  aerosol  mass  spectrometer  (AMS),  single  par¬ 
ticle  soot  photometer  (SP2),  nephelometer,  LiCor  CO2  an¬ 
alyzer,  a  chemiluminescence  ozone  instrument,  and  a  wing- 
mounted  meteorological  probe.  Our  measurements  included: 
CO2;  CO;  NOx;  NH3;  non-methane  organic  compounds;  or¬ 
ganic  aerosol  (OA);  inorganic  aerosol  (nitrate,  ammonium, 
sulfate,  and  chloride);  aerosol  light  scattering;  refractory 
black  carbon  (rBC);  and  ambient  temperature,  relative  hu¬ 
midity,  barometric  pressure,  and  three-dimensional  wind  ve¬ 
locity.  The  molar  ratio  of  excess  O3  to  excess  CO  in  the 
plume  (AO3/ACO)  increased  from  —5.13  (±1.13)  x  10-3 
to  10.2  (±2.16)  x  10-2  in  ~4.5h  following  smoke  emis¬ 
sion.  Excess  acetic  and  formic  acid  (normalized  to  excess 
CO)  increased  by  factors  of  1.73  ±  0.43  and  7.34  ±  3.03  (re¬ 
spectively)  over  the  same  time  since  emission.  Based  on 
the  rapid  decay  of  C2H4  we  infer  an  in-plume  average  OH 


concentration  of  5.27  (±0.97)  x  106  molec  cm-3,  consistent 
with  previous  studies  showing  elevated  OH  concentrations 
in  biomass  burning  plumes.  Ammonium,  nitrate,  and  sulfate 
all  increased  over  the  course  of  4h.  The  observed  ammo¬ 
nium  increase  was  a  factor  of  3.90  ±2.93  in  about  4h,  but 
accounted  for  just  ~36  %  of  the  gaseous  ammonia  lost  on 
a  molar  basis.  Some  of  the  gas  phase  NH3  loss  may  have 
been  due  to  condensation  on,  or  formation  of,  particles  be¬ 
low  the  AMS  detection  range.  NOx  was  converted  to  PAN 
and  particle  nitrate  with  PAN  production  being  about  two 
times  greater  than  production  of  observable  nitrate  in  the  first 
~4  h  following  emission.  The  excess  aerosol  light  scattering 
in  the  plume  (normalized  to  excess  CO2)  increased  by  a  fac¬ 
tor  of  2.50  ±0.74  over  4h.  The  increase  in  light  scattering 
was  similar  to  that  observed  in  an  earlier  study  of  a  biomass 
burning  plume  in  Mexico  where  significant  secondary  forma¬ 
tion  of  OA  closely  tracked  the  increase  in  scattering.  In  the 
California  plume,  however,  AOA/ACO2  decreased  sharply 
for  the  first  hour  and  then  increased  slowly  with  a  net  de¬ 
crease  of  ~20  %  over  4  h.  The  fraction  of  thickly  coated  rBC 
particles  increased  up  to  ~85  %  over  the  4  h  aging  period. 
Decreasing  OA  accompanied  by  increased  scattering/particle 
coating  in  initial  aging  may  be  due  to  a  combination  of  par¬ 
ticle  coagulation  and  evaporation  processes.  Recondensation 
of  species  initially  evaporated  from  the  particles  may  have 
contributed  to  the  subsequent  slow  rise  in  OA.  We  compare 
our  results  to  observations  from  other  plume  aging  studies 
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were  carried  out  onboard  a  Twin  Otter  aircraft  outfitted  with  an  airborne  Fourier  transform  infrared 
spectrometer  (AFTIR),  aerosol  mass  spectrometer  (AMS),  single  particle  soot  photometer  (SP2), 
nephelometer,  LiCor  C02  analyzer  a  chemiluminescence  ozone  instrument,  and  a  wingmounted 
meteorological  probe.  Our  measurements  included  C02;  CO;  NOx;  NH3;  non-methane  organic 
compounds;  organic  aerosol  (OA);  inorganic  aerosol  (nitrate,  ammonium  sulfate,  and  chloride);  aerosol 
light  scattering;  refractory  black  carbon  (rBC);  and  ambient  temperature,  relative  humidity  barometric 
pressure,  and  three-dimensional  wind  velocity.  The  molar  ratio  of  excess  03  to  excess  CO  in  the  plume 
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smoke  emission.  Excess  acetic  and  formic  acid  (normalized  to  excess  CO)  increased  by  factors  of  1.7370.43 
and  7.3473.03  (respectively)  over  the  same  time  since  emission.  Based  on  the  rapid  decay  of  C2H4  we  infer 
an  in-plume  average  OH  concentration  of  5.27  (70.97)7106  molec  cm&#8722;3,  consistent  with  previous 
studies  showing  elevated  OH  concentrations  in  biomass  burning  plumes.  Ammonium,  nitrate,  and  sulfate 
all  increased  over  the  course  of  4  h.  The  observed  ammonium  increase  was  a  factor  of  3.9072.93  in  about  4 
h,  but  accounted  for  just  36%  of  the  gaseous  ammonia  lost  on  a  molar  basis.  Some  of  the  gas  phase  NH3 
loss  may  have  been  due  to  condensation  on,  or  formation  of,  particles  below  the  AMS  detection  range.  NOx 
was  converted  to  PAN  and  particle  nitrate  with  PAN  production  being  about  two  times  greater  than 
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and  suggest  that  differences  in  environmental  factors  such  as 
smoke  concentration,  oxidant  concentration,  actinic  flux,  and 
RH  contribute  significantly  to  the  variation  in  plume  evolu¬ 
tion  observations. 


1  Introduction 

Biomass  burning  (BB)  is  the  largest  global  source  of  primary, 
fine  carbonaceous  particles;  the  second  largest  source  of  trace 
gases;  and  a  key  element  of  the  carbon  cycle  (Crutzen  and 
Andreae,  1990;  Fearnside  et  al.,  1993;  Bond  et  ah,  2004; 
Akagi  et  ah,  2011).  Globally,  dry  season  deforestation  and 
savanna  fires  in  the  tropics  are  the  most  dominant,  and  most 
studied,  type  of  open  burning  (Yokelson  et  ah,  2003a;  van 
der  Werf  et  ah,  2010).  However,  the  comparatively  smaller 
amount  of  BB  that  occurs  year-round  in  the  contiguous  US 
is  a  prominent  domestic  source  of  gas  and  particle  emissions. 
These  fires  may  be  prescribed  or  ignited  naturally.  Prescribed 
burning  is  a  significant  fraction  of  US  BB  and  it  is  used  to  re¬ 
store  or  maintain  the  natural,  beneficial  role  of  fire  in  many 
terrestrial  ecosystems;  reduce  fire  risk  by  consuming  accu¬ 
mulated  wildland  fuels  under  preferred  weather  conditions; 
and  accomplish  other  land  management  objectives  (Biswell, 
1989;  Hardy  et  ah,  2001;  Carter  and  Foster,  2004;  Wied- 
inmyer  and  Hurteau,  2010).  These  burns  are  often  recur¬ 
rent,  sometimes  near  populated  areas,  and  despite  being  con¬ 
trolled  may  occasionally  have  negative  air  quality  impacts. 
Together,  wild  and  prescribed  fires  are  thought  to  produce 
about  one-third  of  the  total  PM2.5  (particle  mass  with  aero¬ 
dynamic  diameter  <2.5  pm)  in  the  US  (Watson,  2002;  Park 
et  ah,  2007).  The  particles  can  impact  local  to  regional  air 
quality,  health,  and  visibility  (McMeeking  et  ah,  2006).  In 
addition,  fires  can  influence  regional  levels  of  ozone  (O3), 
which  is  an  air  toxic  (Pfister  et  ah,  2006).  Large  areas  of 
the  US  commonly  have  O3  mixing  ratios  not  far  below  the 
US  national  ambient  O3  air  quality  standard  (75ppb  8-h  av¬ 
erage).  Thus,  in  these  areas  even  modest  production  of  O3  in 
a  smoke  plume  can  potentially  lead  to  O3  levels  that  exceed 
air  quality  criteria.  Because  fire  has  both  benefits  and  draw¬ 
backs,  the  optimum  amount  and  timing  of  prescribed  fire  is 
an  area  of  active  research  and  better  information  on  the  ini¬ 
tial  chemistry  and  evolution  of  BB  emissions  is  a  major  need 
(Haines  and  Cleaves,  1999;  Sandberg  et  ah,  2002;  Yoder  et 
ah,  2003;  Kauffman,  2004;  Stephens  et  ah,  2007;  Akagi  et 
ah,  2011). 

Next  we  summarize  the  limited  amount  of  data  available 
on  the  emissions  and  downwind  chemical  evolution  of  trace 
gases  and  particles  from  BB  in  the  contiguous  US.  Ground- 
based  measurements  of  prescribed  fire  emissions  from  vari¬ 
ous  temperate  fuel  types  in  the  US  were  carried  out  by  Hardy 
and  Teesdale  (1991).  Hardy  et  ah  (1996)  measured  smoke 
emissions  from  chaparral  fires  in  southern  California  (CA) 
and  reported  emission  factors  (EF,  g  compound  emitted  per 
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kg  dry  biomass  burned)  for  PM  (total  PM,  PM2.5,  PM10), 
CO,  CO2,  CH4,  and  non-methane  hydrocarbons  by  combus¬ 
tion  process  (i.e.  flaming  or  smoldering).  Prescribed  fire 
emissions  have  also  been  measured  from  airborne  platforms 
(Cofer  III.  et  ah,  1988;  Radke  et  ah,  1990;  Yokelson  et  ah, 
1999).  Ward  et  ah  (1993)  compiled  the  first  extensive  in¬ 
ventory  of  particulate  matter  and  trace  gas  EF  for  prescribed 
fires  in  the  US,  which  estimated  EF  based  on  combustion  effi¬ 
ciency  from  field  and  laboratory  measurements  for  the  differ¬ 
ent  classes  of  biomass  consumed.  Numerous  laboratory  stud¬ 
ies  of  the  initial  emissions  from  burning  temperate  ecosystem 
fuels  have  been  conducted  (e.g.  Yokelson  et  ah,  1996,  1997; 
Goode  et  ah,  1999;  McMeeking  et  ah,  2009;  Veres  et  ah, 
2010;  Burling  et  ah,  2010;  Hosseini  et  ah,  2010).  The  stud¬ 
ies  discussed  above  focused  on  measuring  fresh  smoke  emis¬ 
sions,  which  can  help  properly  initialize  models  that  predict 
the  evolution  of  BB  smoke  as  it  ages  downwind.  However, 
they  did  not  directly  observe  post-emission  photochemical 
and  heterogeneous  processes.  Recently,  several  laboratory 
studies  investigated  the  post  emission  processing  of  smoke 
from  temperate  fuels  (Grieshop  et  ah,  2009;  Heringa  et  ah, 
2011;  Hennigan  et  ah,  2011;  Cubison  et  ah,  2011),  but  the  ag¬ 
ing  of  smoke  from  real  fires  is  difficult  to  reproduce  in  a  lab¬ 
oratory.  The  behavior  and  fate  of  these  emissions  are  impor¬ 
tant,  as  air  pollutants  from  fires  can  remain  in  the  atmosphere 
for  many  days  and  they  can  also  be  transported  thousands  of 
kilometers  (Lee  et  ah,  2005;  Yokelson  et  ah,  2007a;  Warneke 
et  ah,  2009).  The  PM2.5  emissions  interact  with  trace  gases 
via  heterogeneous  chemistry  and  gas-particle  interconversion 
processes  as  smoke  plumes  age  (Miihle  et  ah,  2007;  Phule- 
ria  et  ah,  2005).  The  gas-particle  interactions  can  increase 
or  decrease  the  amount  of  PM2.5  and  also  change  the  aerosol 
mixing  state  and  lifetime  (Petters  et  ah,  2009).  Some  trace 
gases  in  BB  smoke  are  rapidly  depleted  post-emission  via  re¬ 
action  with  the  hydroxyl  radical  (OH),  while  others  (e.g.  O3, 
HCHO)  are  created  from  oxidation  of  non-methane  organic 
compounds  (NMOC).  Most  previous  field  studies  of  smoke 
aging  focused  on  major  biomass  burning  regions  such  as  the 
tropics  (Reid  et  ah,  1998;  Hobbs  et  ah,  2003;  Yokelson  et 
ah,  2009),  where  ~70-80  %  of  biomass  burning  is  thought 
to  occur  (Crutzen  and  Andreae,  1990;  van  der  Werf  et  ah, 
2010).  The  only  field  study  that  tracked  the  evolution  of  an 
isolated  aging  plume  from  a  fire  in  temperate  fuels  that  we 
know  of  followed  the  evolution  of  the  plume  from  the  Quinalt 
prescribed  fire  in  the  US  Pacific  Northwest  for  only  30  min 
(Trentmann  et  ah,  2003).  Thus,  more  field  measurements  of 
plume  evolution  in  temperate  ecosystems  are  needed  to  im¬ 
prove  our  understanding  of  the  trade-offs  between  air  quality 
and  ecosystem  health  and  between  wild  and  prescribed  fires 
(Sandberg  et  ah,  2002). 

Two  plume  evolution  processes  are  especially  important, 
both  domestically  and  globally,  but  are  poorly  understood: 
(1)  secondary  organic  aerosol  (SOA)  formation,  and  (2)  O3 
formation.  The  atmospheric  influence  of  pollution  sources  is 
studied  by  modeling  the  dispersion  and  chemical  evolution 
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of  the  emissions  with  local  to  global  scale  models  (e.g.  Al¬ 
varado  and  Prinn,  2009;  Byim  and  Schere,  2006;  Heald  et  al., 
2005).  Recently,  numerous  studies  in  urban  plumes  have  ob¬ 
served  rapid  production  of  SOA  in  amounts  up  to  ten  times 
the  original  emissions  of  primary  organic  aerosol  (POA), 
which  greatly  exceeds  the  predictions  of  standard  aerosol 
modules  in  atmospheric  chemistry  models  (de  Gouw  and 
Jimenez,  2009).  This  has  spurred  the  development  of  new, 
largely  empirical  modeling  approaches  that  can  cope  with  the 
large  increases  (Robinson  et  al.,  2007).  In  contrast,  measure¬ 
ments  of  organic  aerosol  (OA)  evolution  in  biomass  burning 
smoke  are  rare  and  have  thus  far  shown  smaller  and  variable 
SOA  formation  compared  to  urban  plumes.  In  an  airborne 
study  of  a  BB  plume  in  Mexico  (Yokelson  et  al.,  2009)  and 
a  laboratory  study  of  wood  smoke  (Grieshop  et  al.,  2009) 
the  OA/POA  ratio  increased  to  greater  than  two  in  about  1 
and  3  h,  respectively.  Capes  et  al.  (2009)  did  not  observe  evi¬ 
dence  of  SOA  formation  in  African  BB  emissions  on  the  time 
scale  of  several  days.  Cubison  et  al.  (201 1)  observed  modest 
SOA  formation  that  averaged  about  25  %  of  POA  in  boreal 
BB  smoke.  The  laboratory  study  of  Hennigan  et  al.  (2011) 
observed  OA/POA  between  0.7  and  2.9  in  BB  smoke  aged 
for  3 — 4.5  h,  while  Heringa  et  al.  (2011)  found  that  OA/POA 
rose  to  4.1  ±  1.4  after  5  h  for  wood  smoke  aged  in  a  smog 
chamber.  We  also  note  that  biomass  burning  is  known  to 
emit  species  that  are  precursors  for  aqueous  phase  SOA  such 
as  glycolaldehyde  and  phenol  (Yokelson  et  al.,  1997;  Perri  et 
al.,  2009;  Lim  et  al.,  2010)  and  that  a  large  fraction  of  the 
organic  gases  emitted  by  BB  are  still  unidentified  (Karl  et 
al.,  2007;  Warneke  et  al.,  2011)  and  may  play  a  role  in  sev¬ 
eral  types  of  SOA  formation.  Many  more  measurements  of 
OA  evolution  in  BB  plumes  are  needed  to  determine  the  rela¬ 
tive  importance  of  the  different  OA/POA  outcomes  at  various 
time  scales  and  to  understand  the  factors  driving  the  variabil¬ 
ity. 

A  second  major  issue  is  understanding  O3  production  in 
biomass  burning  plumes.  O3  production  is  evidently  the 
norm  in  tropical  BB  plumes,  but  both  production  and  destruc¬ 
tion  of  O3  occurs  in  boreal  BB  plumes  and  little  is  known 
about  O3  chemistry  in  temperate  BB  plumes  (Akagi  et  al., 
2011).  The  few  available  time-resolved  observations  of  O3 
formation  in  isolated  plumes  have  proven  difficult  to  repro¬ 
duce  even  in  specialized  photochemical  models  (Trentmann 
et  al.,  2005;  Alvarado  and  Prinn,  2009).  Additional  mea¬ 
surements  of  the  O3  formation  rate  in  various  plume  environ¬ 
ments  are  needed  to  further  develop  the  specialized  models 
and  ultimately  guide  the  implementation  of  reasonably  accu¬ 
rate  O3  formation  mechanisms  into  air  quality  models  cur¬ 
rently  in  wide  use,  such  as  CMAQ  (Byun  and  Schere,  2006). 

Here  we  report  new  field  measurements  that  address  the  is¬ 
sues  discussed  above.  We  acquired  the  measurements  during 
the  San  Luis  Obispo  Biomass  Burning  (SLOBB)  experiment, 
which  was  carried  out  in  November  2009  primarily  in  mid¬ 
coast  California.  The  main  goal  of  SLOBB  was  to  study  the 
emissions  and  the  evolution  of  gases  and  particles  produced 
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by  prescribed  burns  on  US  Department  of  Defense  lands.  We 
measured  many  new  emission  factors  for  US  prescribed  fires 
during  SLOBB  and  related  field  studies  as  detailed  by  Burl¬ 
ing  et  al.  (2011).  In  addition,  a  prescribed  fire  near  Buell- 
ton,  CA  named  the  Williams  Fire  proved  particularly  oppor¬ 
tune  for  the  study  of  gas-phase  photochemistry  and  particle 
transformations  in  a  smoke  plume.  These  processes  are  dif¬ 
ficult  to  simulate  in  a  laboratory  and  frequently  inaccessible 
even  for  airborne  platforms.  The  Williams  Fire  produced  a 
long,  clearly-identifiable,  isolated  plume  above  the  boundary 
layer  where  the  smoke  evolved  without  influence  from  fac¬ 
tors  such  as  clouds,  mixing  with  plumes  from  other  pollution 
sources,  or  complex  meteorology.  Thus,  the  measurements 
provide  an  ideal  case  study  to  test  the  chemical  mechanism 
of  a  smoke  model.  To  our  knowledge  the  Williams  Fire  is  the 
first  temperate  zone  fire  for  which  the  evolution  of  numerous 
trace  gases  and  aerosol  species  in  the  plume  was  continu¬ 
ously  tracked  for  smoke  ages  of  more  than  one -half  hour.  Ex¬ 
tensive  particle  and  trace  gas  measurements  in  fresh  and  aged 
smoke  (for  ages  up  to  4.5  h  old)  were  performed  to  quantify 
the  changes  in  NMOC,  O3,  NOx,  PAN,  NH3,  light  scatter¬ 
ing,  inorganic  and  organic  aerosol,  refractory  black  carbon 
(rBC),  and  the  rBC  mixing  state.  We  describe  and  discuss 
these  measurements  in  detail  in  the  following  sections. 

2  Experimental  details 

2.1  Instrument  details 

2.1.1  Airborne  Fourier  transform  infrared 
spectrometer  (AFTIR) 

The  University  of  Montana  AFTIR  system  and  the  instru¬ 
ments  described  below  were  deployed  on  a  US  Forest  Ser¬ 
vice  (USFS)  Twin  Otter.  The  AFTIR  was  used  to  measure 
21  gas-phase  species:  water  vapor  (H2O),  carbon  dioxide 
(CO2),  carbon  monoxide  (CO),  methane  (CH4),  nitric  oxide 
(NO),  nitrogen  dioxide  (NO2),  ammonia  (NH3),  hydrogen 
cyanide  (HCN),  nitrous  acid  (HONO),  peroxy  acetyl  nitrate 
(PAN),  ozone  (O3),  glycolaldehyde  (HCOCH2OH),  ethy¬ 
lene  (C2H4),  acetylene  (C2H2),  propylene  (C3H6),  formalde¬ 
hyde  (HCHO),  methanol  (CH3OH),  furan  (C4H4O),  phe¬ 
nol  (C6H5OH),  acetic  acid  (CH3COOH),  and  formic  acid 
(HCOOH).  Ram  air  was  directed  through  a  halocarbon-wax 
coated  inlet  and  into  a  Pyrex  multipass  cell.  The  IR  beam  was 
directed  into  the  cell  where  it  traversed  a  total  path  length  of 
78  m  and  was  then  focused  onto  an  MCT  detector.  The  cell 
exchange  time  was  about  ten  seconds  when  the  flow  con¬ 
trol  valves  were  open  and  IR  spectra  were  collected  at  1  Hz. 
“Grab  samples”  of  air  were  selected  by  closing  the  valves  for 
1-2  min  to  temporarily  stop  system  flow  to  allow  signal  aver¬ 
aging  on  the  smoke  or  background  samples.  The  IR  spectra 
were  analyzed  to  identify  and  quantify  all  detectable  com¬ 
pounds  (Sect.  2.4.1).  More  details  of  the  AFTIR  system  are 
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given  by  Yokelson  et  al.  (1999,  2003a).  In  this  project,  up¬ 
grades  to  the  spectrometer  electronics,  transfer  optics,  and 
multipass  cell  led  to  a  factor  of  three  increase  in  sensitivity 
compared  to  the  original  system. 

2.1.2  Aerosol  mass  spectrometer  (AMS) 

An  Aerodyne  compact  time-of-flight  (CToF)  aerosol  mass 
spectrometer  (herein  referred  to  as  AMS)  measured  aerosol 
chemical  composition  in  a  repeating  cycle  for  4  out  of  ev¬ 
ery  12  s  during  flight,  including  within  the  smoke  plume.  An 
isokinetic  particle  inlet  sampling  fine  particles  with  a  diam¬ 
eter  cut-off  of  a  few  microns  (Yokelson  et  al.,  2007b)  sup¬ 
plied  the  AMS  and  a  single  particle  soot  photometer  (SP2, 
Sect.  2.1.3).  The  AMS  and  SP2  do  not  measure  super-micron 
particles,  so  the  inlet  transmission  should  not  have  affected 
our  results.  In  addition,  particles  smaller  than  1  pm  diame¬ 
ter  account  for  nearly  all  the  fine  particle  mass  emitted  by 
biomass  fires  (Radke  et  al.,  1991;  Reid  et  al.,  2005b),  so  our 
composition  analyses  for  fine  particles  should  not  be  affected 
by  the  lack  of  sensitivity  to  super-micron  particles.  The  AMS 
collected  sub-micron  particles  via  an  aerodynamic  lens  into 
a  high  vacuum  particle  sizing  chamber.  At  the  end  of  the 
particle  sizing  chamber,  the  particles  impact  a  600  °C  vapor¬ 
izer  and  filament  assembly  where  they  are  vaporized  and  ion¬ 
ized  by  electron  impact.  The  resulting  molecular  fragments 
are  then  extracted  into  an  ion  time-of-flight  chamber  where 
they  are  detected  and  interpreted  as  mass  spectra.  The  AMS 
has  been  described  in  great  detail  elsewhere  (Drewnick  et  al., 
2005;  Canagaratna  et  al.,  2007).  A  collection  efficiency  of 
0.5  (Huffman  et  al.,  2005;  Drewnick  et  al.,  2003;  Allan  et  al., 
2004)  was  applied  to  the  AMS  data,  which  were  processed 
with  the  ToF-AMS  Analysis  Toolkit  in  Igor  Pro  6  (Wavemet- 
rics.  Lake  Oswego,  OR)  to  retrieve  the  mass  concentration  at 
standard  temperature  and  pressure  (pgsm-3,  273  K,  1  atm) 
for  the  major  non-refractory  particle  species;  OA,  non-sea 
salt  chloride,  nitrate,  sulfate,  and  ammonium,  with  <36% 
uncertainty.  Adjustments  to  the  default  AMS  unit  mass  reso¬ 
lution  (UMR)  fragmentation  table  were  made  for  sulfate  and 
nitrate  concentrations.  For  instance,  nitrate  was  corrected  for 
organic  interference  at  the  nitrate  peak  at  m/z  30  based  on  the 
results  of  Bae  et  al.  (2007). 

2.1.3  Single  particle  soot  photometer  (SP2) 

Refractory  black  carbon  (rBC,  pgsm  3)  was  measured  us¬ 
ing  a  single  particle  soot  photometer  (SP2)  (Stephens  et  al., 
2003).  The  isokinetic  particle  inlet  that  supplied  the  AMS 
was  followed  by  a  three-way  splitter  and  also  supplied  the 
SP2  with  carefully  matched  flow  rates.  Particles  were  then 
drawn  through  a  1064  nm  Nd:YAG  laser  cavity  where  light 
scattered  by  the  particles  was  measured  by  two  avalanche 
photodiode  detectors.  Sufficiently  light-absorbing  particles 
(at  1064  nm)  were  heated  and  eventually  vaporized  when 
passing  through  the  beam.  The  emitted  incandescent  light 
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was  detected  by  two  photomultiplier  tubes  with  optical  fil¬ 
ters  to  measure  light  over  two  wavelength  ranges.  The  in¬ 
candescence  signal  was  related  linearly  to  the  mass  of  rBC, 
regardless  of  the  presence  or  type  of  rBC  coating  (Moteki 
and  Rondo,  2007).  The  sample  flow  into  the  SP2  was  diluted 
with  filtered  air  at  a  ratio  of  7:1  for  smoke  plume  samples 
near  the  source  with  high  concentrations  of  particles.  A  di¬ 
lution  of  0.35:1  or  no  dilution  was  used  for  the  less  concen¬ 
trated  downwind  samples.  Information  on  the  rBC  mixing 
state  was  inferred  from  the  time  delay  between  the  peak  of 
the  scattering  and  incandescence  signals,  which  is  a  proxy 
for  the  time  taken  for  coatings  to  evaporate  before  incandes¬ 
cence  can  occur  (Schwarz  et  al.,  2006).  Particles  with  long 
time  delays  were  considered  “thickly  coated”,  while  those 
with  shorter  delays  were  “thinly  coated"  following  Shiraiwa 
et  al.  (2007).  The  critical  delay  time  needed  for  classifica¬ 
tion  was  based  on  the  observed  minimum  in  the  bimodal  fre¬ 
quency  distribution  of  delay  times  for  the  rBC  particles  be¬ 
tween  120  and  315  nm  volume  equivalent  diameter  and  was 
1.1  ps.  The  SP2  data  products  were  collected  at  >1  Hz  and 
later  averaged  as  needed  to  match  the  sampling  rates  of  the 
other  instruments. 

2.1.4  Other  simultaneous  measurements 

A  25mmi.d.  forward  facing  elbow  “fastflow”  inlet,  co¬ 
located  with  the  isokinetic  and  AFTIR  inlets,  fed  air  to 
a  Radiance  Research  Model  903  nephelometer  fit  with  a 
2.5  pm  cut-off  cyclone  that  measured  the  dry  particle  light 
scattering  coefficient  (/?scat)  at  530  nm  at  0.5  Hz.  The  /;scat 
was  converted  to  £>scat  at  STP  and  then  to  a  mass  of  PM2.5 
(pgsm-3)  using  a  conversion  factor  (for  fresh  smoke)  of 
208  800  pg  sm-2  (equivalent  to  a  mass  scattering  efficiency 
of  ~4.8sm2g-1)  obtained  by  a  gravimetric  calibration  of 
the  type  described  by  Trent  et  al.  (2000).  The  nephelome¬ 
ter  inlet  also  supplied  air  to  a  non-dispersive  infrared  instru¬ 
ment  (NDIR,  LiCor  model  6262),  which  measured  CO2  at 
the  same  time  as  £>SCat-  Since  the  AFTIR  also  measured  CO2, 
this  allowed  the  PM2.5  measurements  to  be  normalized  to 
the  trace  gas  data  as  described  in  detail  below.  A  second 
NDIR  (LiCor  model  7000)  measured  CO2  (at  0.5  Hz)  from 
the  third  channel  on  the  isokinetic  particle  inlet  that  also  sup¬ 
plied  the  AMS  and  SP2  allowing  the  data  from  the  latter  in¬ 
struments  to  be  coupled  to  the  trace  gas  data.  Measurements 
of  the  ambient  three-dimensional  wind  velocity,  tempera¬ 
ture,  relative  humidity,  and  barometric  pressure  at  1-Hz  were 
obtained  with  a  wing-mounted  Aircraft  Integrated  Meteoro¬ 
logical  Measuring  System  probe  (AIMMS-20,  Aventech  Re¬ 
search,  Inc.).  Details  of  the  AIMMS-20  probe  and  a  perfor¬ 
mance  evaluation  are  provided  by  Beswick  et  al.  (2008).  A 
chemiluminescence  instrument  measured  ozone  at  0.5  Hz. 
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Fig.  1.  (a)  Photograph  of  the  high-  and  low-altitude  smoke  plumes  from  the  Williams  Fire  taken  at  12:20  LT  on  17  November  2009  from 
the  USFS  Twin  Otter,  (b)  HYSPLIT  forward  trajectories  starting  at  11:00LT  at  the  Williams  Fire  (Draxler  and  Rolph,  2010;  Rolph,  2010). 
Eight  24  h  trajectories  show  air  at  600-1000  ma.m.s.l.  (green)  drifting  SE  and  air  at  1 100-1700 ma.m.s.l.  (purple)  drifting  NE.  (c)  Vertical 
profile  of  windspeeds  from  Flights  1  (green)  and  2  (blue). 


2.2  Vegetation  and  fire  development 

The  Williams  prescribed  fire  was  carried  out  on  17  November 
2009.  The  81  hectare  site  (~10km  north  of  Buellton,  Cal¬ 
ifornia;  34°4T45"N,  120°12'23"W)  covered  several  hill¬ 
sides  with  various  aspects.  The  land  cover  can  be  classified 
as  chaparral  and  the  vegetation  burned  consisted  of  scrub 
oak  woodland  understory  and  coastal  sage  scrub,  which  is 
a  mixture  of  California  sagebrush  ( Artemisia  californica) 
with  coyote  brush  ( Baccharis  pilularis ),  black  sage  ( Salvia 
mellifera),  and  California  buckwheat  (Eriogonum  fascicula- 
turri).  Dead  fuel  moisture  ranged  from  9  to  15  %  with  higher 
moisture  in  the  larger  diameter  fuels.  Live  fuel  moisture 
was  70%,  which  is  below  average  compared  to  normally 
observed  values.  Ground  level  winds  were  light  and  vari¬ 
able  at  2-8  km  h-1.  Ground  level  relative  humidity  (RH) 
ranged  from  26-11  %  between  09:00-1 3:00  LT  (PST,  UTC- 
8).  Surface  temperature  ranged  from  19  °C  at  09:00  LT  to 
24  °C  at  12:20  LT.  Skies  were  clear  all  day.  (Fuels  and  sur¬ 
face  weather  information,  J.  McCarty,  personal  communica¬ 
tion,  Santa  Barbara  County  Fire  Department,  2009.)  Igni¬ 


tion  of  the  Williams  Fire  began  at  approximately  09:57  LT. 
The  plume  built  up  gradually  and  most  of  the  smoke  rose 
to  ~  1200-1 336  m  above  mean  sea  level  (ma.m.s.l.)  and 
then  drifted  northeast.  Some  smoke  detrained  at  lower  al¬ 
titudes  (below  ~600  ma.m.s.l.)  and  drifted  southeast,  set¬ 
tling  into  nearby  valleys  (Fig.  la).  Our  visual  observa¬ 
tions  were  consistent  with  air  parcel  forward  trajectory  sim¬ 
ulations  using  HYSPLIT  (NOAA  Air  Resources  Labora¬ 
tory,  http://ready.arl.noaa.gov/HYSPLIT.php),  which  shows 
low-altitude  trajectories  (600-1000 ma.m.s.l..  Fig.  lb)  orig¬ 
inating  at  the  source  traveling  southeast,  while  higher  al¬ 
titude  trajectories  (1100-1700ma.m.s.l.)  traveled  north¬ 
east  (Fig.  lb).  The  upwind  vertical  profiles  of  temperature 
and  water  vapor  mixing  ratio  from  the  AIMMS-20  (Sup¬ 
plement  Fig.  SI)  showed  that  the  top  of  the  boundary  layer 
was  located  at  ~800 ma.m.s.l.  (~450m  above  ground  level, 
or  a.g.l.)  and  the  AIMMS-20  wind  direction  data  con¬ 
firmed  the  shift  in  wind  direction  modeled  by  HYSPLIT 
at  1 100 i  100 ma.m.s.l.  The  bulk  of  the  smoke  was  suffi¬ 
ciently  buoyant  to  penetrate  the  temperature  inversion  and 
rose  to  ~1200  ma.m.s.l.  where  it  encountered  an  average 
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Fig.  2.  (a)  Side  view  of  the  Williams  Fire  high  altitude  plume  at  12:40  LT  from  ~58km  downwind,  (b)  HYSPLIT  trajectory  for  source  air 
starting  at  11:00  LT  and  1200ma.m.s.l.  No  significant  altitude  gain  is  modeled  in  the  first  ~4hsince  emission. 


windspeed  of  ~5ms_1  (Fig.  lc).  Low-altitude  smoke  dis¬ 
persed  slowly  (~2.5  ms-1  was  the  average  windspeed  at  al¬ 
titudes  <1200ma.m.s.l.  during  Flights  1  and  2)  and  tended 
to  collect  in  low-lying  areas  as  discussed  further  in  Sect.  3.5. 

The  high  altitude  plume  that  drifted  northeast  from  the 
source  was  isolated  from  other  pollution  sources  and  unim¬ 
peded  by  topography  as  seen  in  the  photograph  taken 
from  ~58  km  downwind  (Fig.  2a).  The  high  altitude 
plume  gradually  rose  from  ~1200ma.m.s.l.  at  the  source 
to  ~1 800  m  a.m.s.l.  at  ~58km  (4.5  h)  downwind,  but  HYS¬ 
PLIT  trajectories  starting  at  1200ma.m.s.l.  failed  to  cap¬ 
ture  this  observed  behavior,  showing  no  substantial  altitude 
gain/loss  for  at  least  4  h  following  fire  ignition  (Fig.  2b).  The 
observed  rise  might  have  been  due  to  self  lofting  since  the 
high  levels  of  O3  and  black  carbon  in  the  plume  would  read¬ 
ily  absorb  significant  amounts  of  solar  energy.  By  14:30 LT 
the  fire  had  transitioned  to  mainly  smoldering  combustion 
with  minimal  plume  rise. 

2.3  Flight  details 

Two  flights  (mid-morning  and  afternoon)  were  conducted  on¬ 
board  the  Twin  Otter  aircraft  to  sample  initial  emissions  and 
aged  smoke.  Table  1  summarizes  the  times,  locations,  and 
meteorological  conditions  for  AFTIR  samples  taken  during 
both  flights.  For  every  plume  sample  on  both  flights  we  also 
collected  a  background  sample  at  similar  altitude  just  outside 
the  plume. 

Flight  1  occurred  from  10:29-13:10LT.  We  arrived  at  the 
site  of  the  fire  approximately  30  min  after  initial  ignition,  at 
which  time  the  fire  was  beginning  to  produce  a  large  plume 
of  thick,  white  and  brown  smoke.  From  10:50-12: 19  LT  we 
collected  ten  samples  of  “nascent”  smoke  (smoke  “just  cre¬ 
ated”  or  “emerging  from  another  form”)  (Samples  1-6,  8, 
10-12  in  Table  1).  The  source  samples  were  collected  ~600- 
1300  m  a.m.s.l.  at  the  center  and  top  portions  of  the  single  up¬ 
draft  core.  After  12:30  LT  we  concentrated  on  sampling  the 
smoke  downwind  at  ~1200m  a.m.s.l.,  making  perpendicular 


plume  transects  in  an  “S”  shaped  pattern  and  sampling  the 
smoke  eight  times  4.5-16.  lkm  from  the  source  (see  flight 
path  and  AFTIR  smoke  sample  locations  in  Fig.  3). 

Flight  2  was  conducted  from  14:20-15:47  LT  and  sam¬ 
pled  the  high  altitude  plume  at  1200-1900  m  a.m.s.l.  and 
27-58  km  from  the  source.  From  14:45-15: 16  LT  we  col¬ 
lected  10  grab  samples  of  aged  smoke  2. 4-4. 5  h  old.  We  ob¬ 
tained  estimated  emission  times  by  first  calculating  the  aver¬ 
age  windspeed  for  incremental  altitude  bins  of  100  m  a.m.s.l. 
The  smoke  sample  distance  from  the  plume  source  was  then 
divided  by  the  average  windspeed  at  the  sample  altitude.  Es¬ 
timated  times  since  emission,  or  smoke  “ages”,  were  calcu¬ 
lated  for  all  the  downwind  samples  from  both  flights  (Ta¬ 
ble  1)  using  this  methodology.  This  method  implies  that  the 
oldest  smoke  sample  had  aged  approximately  4.5  h.  In  ad¬ 
dition,  based  on  the  estimated  emission  time  of  each  down¬ 
wind  sample:  14  of  the  18  downwind  samples  from  the  two 
flights  were  likely  emitted  during  the  time  we  were  sampling 
the  source  (10:50-12: 19  LT)  and  the  other  4  downwind  sam¬ 
ples  were  emitted  within  20  min  of  our  source  sampling  pe¬ 
riod,  supporting  that  the  data  are  pseudo-Lagrangian.  The 
Williams  Fire  plume  was  visible  in  GOES  images  and  the 
position  of  the  plume  in  the  last  daytime  image  (~  16:45  LT, 
acquired  after  Flight  2)  is  indicated  in  yellow  in  Fig.  3.  Our 
measurements  covered  ~70  %  of  the  distance  the  plume  trav¬ 
eled  the  first  day. 

2.4  Data  processing  and  synthesis 
2.4.1  FTIR  spectral  analysis 

We  use  nonlinear  least  squares,  multi-component  fits  to  re¬ 
gions  of  the  IR  spectra  with  the  HITRAN  spectral  database 
(Rothman  et  al.,  2009)  and  reference  spectra  from  the  Pacific 
Northwest  National  Laboratory  (Sharpe  et  al.,  2004;  John¬ 
son  et  al.,  2006,  2010)  to  obtain  mixing  ratios  in  both  the 
smoke  and  background  spectra  for  most  compounds  (Grif¬ 
fith,  1996;  Burling  et  al.,  2010).  For  these  compounds  the 
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Table  1.  Smoke  samples  from  the  Williams  Fire  obtained  during  Flights  1  and  2. 


Sample 

Time 

(LT) 

Latitude 
(°  N) 

Longitude 

(°W) 

Altitude 

(ma.m.s.l.) 

Pressure 

(hPa) 

Relative 
Humidity 
(RH  %) 

Outside  Air 
Temperature 
(°C) 

Distance 
From  Fire 
Source  (km) 

Time  Since 
Emission 
(h:min:s) 

Estimated 
Emission 
Time  (LT) 

Flight  1 

1 

10:51:13 

34°  41 '40" 

120°  12' 13" 

817 

923.0 

13.0 

16.6 

0.31 

0:02:04 

10:49:09 

2 

10:55:45 

34°41'41" 

120°  12' 18" 

819 

919.2 

14.0 

17.3 

0.21 

0:01:25 

10:54:20 

3 

11:02:01 

340  41 '42" 

120°  12' 18" 

871 

918.0 

12.4 

17.3 

0.16 

0:01:05 

11:00:56 

4 

11:07:39 

34°  41 '42" 

120°  12'23" 

676 

938.8 

14.2 

17.2 

0.12 

0:00:59 

11:06:40 

5 

11:11:54 

340  41 '42" 

120°12'28" 

748 

925.9 

15.9 

17.9 

0.18 

0:00:58 

11:10:56 

6 

11:20:47 

34°42'4" 

120°  12' 11" 

1336 

865.3 

11.2 

15.4 

0.65 

0:03:05 

11:17:42 

7 

11:31:55 

34°46'35" 

120°  10'59" 

1230 

877.6 

15.5 

15.3 

9.21 

0:37:27 

10:54:28 

8 

11:55:25 

34°42'39" 

120°  12'41" 

1196 

882.3 

13.4 

15.2 

1.71 

0:07:35 

11:47:50 

9 

11:59:58 

34°43'42" 

120°10'35" 

1192 

881.8 

16.8 

15.3 

4.55 

0:20:08 

11:39:50 

10 

12:08:08 

340  41 '44" 

120°  12'  30" 

621 

938.7 

26.3 

17.8 

0.18 

0:01:23 

12:06:45 

11 

12:15:24 

340  41 '45" 

120°  12'23" 

624 

936.3 

25.9 

18.3 

0.00 

0:00:00 

12:15:24 

12 

12:18:55 

34°  42' 6" 

120°  12'33" 

901 

911.4 

23.6 

16.5 

0.68 

0:04:17 

12:14:38 

13 

12:30:30 

34°45'28" 

120°  9' 22" 

1198 

880.3 

15.4 

15.5 

8.28 

0:36:38 

11:53:52 

14 

12:34:01 

34°46'40" 

120° 11 '48" 

1202 

879.9 

15.0 

15.7 

9.15 

0:37:13 

11:56:48 

15 

12:38:50 

34°49'28" 

120°  9' 5  8" 

1195 

879.9 

15.6 

15.4 

14.8 

1:05:26 

11:33:24 

16 

12:44:57 

34°52'18" 

120°  10'44" 

1196 

879.9 

14.4 

15.6 

19.7 

1:27:21 

11:17:36 

17 

12:47:56 

34°50'50" 

120°  6'  26" 

1210 

884.0 

12.5 

15.8 

19.1 

1:17:44 

11:30:12 

18 

12:55:27 

34°48'2l" 

120°10'57" 

1163 

888.6 

16.3 

15.9 

12.4 

0:55:03 

12:00:24 

Flight  2 

19 

14:45:06 

34°55'28" 

120°  6' 15" 

1259 

879.9 

14.7 

15.4 

27.1 

3:17:12 

11:27:54 

20 

14:49:08 

35°4'34" 

120°6'40" 

1744 

825.8 

12.4 

14.0 

43.2 

3:03:30 

11:45:38 

21 

14:50:06 

35°6'3l" 

120°7'33" 

1790 

821.2 

15.3 

12.7 

46.5 

3:17:35 

11:32:31 

22 

14:52:41 

35°8'16" 

120°5'52" 

1898 

810.1 

13.7 

12.9 

50.2 

3:20:17 

11:32:24 

23 

14:56:11 

35°9'40" 

120°  2' 5  6" 

1899 

813.2 

13.5 

12.7 

53.7 

3:34:29 

11:21:42 

24 

14:56:48 

35°10'52" 

120°3'42" 

1882 

813.3 

14.3 

12.2 

55.6 

3:41:56 

11:14:52 

25 

14:59:48 

35°11'56" 

120°2'37" 

1979 

801.2 

13.1 

12.3 

57.9 

4:29:44 

10:30:04 

26 

15:07:59 

35°5'18" 

120°  6' 7" 

1423 

867.7 

14.2 

13.9 

44.7 

2:26:20 

12:41:39 

27 

15:12:00 

34°57'51" 

120°3'46" 

1377 

864.7 

20.3 

13.9 

32.6 

2:43:56 

12:28:04 

28 

15:16:26 

34°57'28" 

120°  4' 13" 

1376 

864.3 

16.6 

13.9 

31.7 

2:39:22 

12:37:04 

excess  mixing  ratio  for  any  species  “A”  (denoted  AA)  is 
the  mixing  ratio  of  species  “A”  in  the  plume  sample  minus 
its  mixing  ratio  in  the  background  air.  For  HCN,  NO,  and 
NO2  we  determine  the  excess  mixing  ratios  directly  from  ab¬ 
sorbance  spectra  via  spectral  subtraction  or  peak  integration 
(Yokelson  et  al.,  1997).  The  absorbance  spectra  are  made  us¬ 
ing  background  spectra  acquired  at  nearly  the  same  pressure. 
The  NO  detection  limit  is  ~25  ppb  and  all  excess  mixing  ra¬ 
tios  below  this  concentration  were  omitted.  The  detection 
limit  for  the  other  compounds  reported  here  is  1-10  ppb. 

2.4.2  Estimation  of  excess  mixing  ratios,  normalized  ex¬ 
cess  mixing  ratios  (NEMR),  emission  ratios  (ER), 
and  emission  factors  (EF) 

The  procedure  to  calculate  excess  mixing  ratios  for  the  AF- 
T1R  data  is  given  above.  Excess  mixing  ratios  can  also 
be  calculated  from  the  continuous  data  (SP2,  CO2,  Z?sc at) 
and  repetitive  discrete  samples  (AMS).  In  both  these  cases. 


we  compute  the  background  value  from  the  average  of  a 
large  number  of  points  at  the  plume  altitude  (but  outside  the 
plume)  and  then  subtract  that  background  from  the  values 
obtained  in  the  plume. 

The  normalized  excess  mixing  ratio  (NEMR)  is  obtained 
by  dividing  AA  by  the  excess  mixing  ratio  of  another  species 
(AT)  measured  at  the  same  time  and  location;  usually  a  fairly 
long-lived  plume  “tracer”  such  as  ACO  or  ACO2.  For  the 
Williams  Fire,  the  NEMRs  for  all  species  (except  ANH3  ref¬ 
erenced  to  ACO  or  ACO2)  measured  at  the  source  of  the 
fire  did  not  vary  greatly  as  the  fire  progressed  through  the 
available  fuels.  We  incorporated  the  ten  measurements  of 
the  source  NEMR  into  an  estimate  of  the  fire-average  ini¬ 
tial  molar  emission  ratio  (ER)  as  described  below.  The  ER 
has  two  important  uses:  (1)  Since  the  CO  or  CO2  tracers  di¬ 
lute  at  the  same  rate  as  the  other  species,  differences  between 
the  ERs  and  the  NEMRs  measured  downwind  allow  us  to 
quantify  post-emission  gas-phase  chemistry  and  gas-particle 
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•  Williams  Fire  Source 

+  Smoke  Samples 

NORTH 

Flight  2  1 420- 1 547  LT 

Williams  Fire  Smoke  1645  LT(GOFS) 

0  5  10  20  30 

40  50 

Kilometers 

Fig.  3.  The  flight  paths  and  the  location  of  smoke  samples  for  Flight  1  (blue)  and  Flight  2  (orange).  The  approximate  position  of  the 
densest  smoke  in  the  Williams  Fire  plume  (yellow)  during  the  last  daytime  GOES  image  (~16:45  LT)  captured  after  our  second  flight  on 
17  November  2009.  The  overlay  shows  that  our  measurements  probed  most  of  the  daytime  plume  evolution  that  occurred  on  17  November 
(plume  image  from  http://www.firedetect.noaa.gov). 


partitioning;  and  (2)  The  ER  can  be  used  to  calculate  emis¬ 
sion  factors  (EF)  as  described  below. 

The  measurement  of  a  NEMR  is  most  accurate  when  both 
species  were  measured  simultaneously  on  the  same  instru¬ 
ment  (e.g.  two  gases  in  the  AFTIR  cell  or  two  particle  species 
on  the  AMS).  When  comparing  NEMR  that  were  determined 
using  data  from  two  different  instruments,  error  can  be  intro¬ 
duced  due  to  the  different  time  responses  of  the  instruments. 
However,  this  error  can  be  largely  eliminated  for  continuous 
instruments  by  deriving  the  NEMR  from  a  comparison  of  the 
integrals  over  the  whole  plume  sample.  In  addition,  in  the 
dispersed,  downwind  plume,  the  excess  mixing  ratios  tended 
to  vary  slowly  in  time  and  space  compared  to  measurement 
frequency.  In  this  work,  we  most  often  use  NEMR  relative 
to  ACO  since  CO  has  a  lifetime  of  several  months  at  ~35°  N 
in  November,  low  background  variability,  and  there  were  no 
other  major  sources  of  CO  noted  in  the  vicinity.  In  some 
cases  we  show  downwind  NEMR  to  CO2  since  that  was  the 
relevant  co-measured  species. 

We  computed  fire-average  initial  emission  ratios  as  our 
best  estimate  of  the  starting  conditions  in  the  plume  using  our 
ten  source  samples  as  follows.  For  grab  samples  we  take  the 
fire-average  ER  as  the  slope  of  the  linear  least-squares  line 
with  the  intercept  forced  to  zero  when  plotting  AX  against 
AT  (Yokelson  et  al.,  1999).  For  example,  the  slope  of  the 
plot  of  AX  vs.  ACO  yields  the  initial  molar  emission  ratio 
for  compound  “A”  to  CO  and  the  uncertainty  is  the  1  a  stan¬ 
dard  error  in  the  slope  (Fig.  4a,  b).  The  intercept  is  forced 


to  zero  since  the  background  concentration  is  typically  well 
known  and  variability  in  the  plume  can  affect  the  intercept  if 
it  is  not  forced.  This  method  heavily  weights  the  large  values 
that  may  reflect  more  intense  combustion  and  were  collected 
with  higher  signal-to-noise.  Smoke  samples  collected  more 
than  1.8  km  from  the  source  showed  signs  of  aging  (O3  for¬ 
mation)  and  were  not  used  to  derive  the  initial  ER.  We  ob¬ 
tained  initial  mass  emission  ratios  for  the  AMS  species,  rBC, 
and  PM2.5  to  CO2  since  CO2  was  measured  on  the  same  in¬ 
let.  These  ratios  were  computed  by  comparing  the  integrated 
excess  amounts  as  the  aircraft  passed  through  the  nascent 
smoke  as  mentioned  above.  For  the  set  of  passes  through  the 
fresh  smoke  we  plot  the  integrals  of  the  excess  values  from 
the  AMS  and  continuous  instruments  versus  the  integrals  of 
excess  CO2  and  obtain  the  ER  from  the  slope,  analogous  to 
the  grab  sample  plots  (Fig.  4c,  d).  To  test  the  robustness  of 
our  integral-based  AMS  ERs  to  CO2,  we  compared  them  to 
ER  obtained  by  plotting  the  excess  values  from  the  4  s  dis¬ 
crete  samples  versus  the  values  of  ACO2  measured  during 
the  same  4  s  interval.  In  nearly  all  cases  we  obtained  agree¬ 
ment  within  5  %,  except  for  sulfate,  which  was  reconstructed 
based  on  samples  of  smoke  with  low  levels  of  organic  load¬ 
ing  where  the  AMS  fragmentation  table  gave  reasonable  re¬ 
sults.  The  major  sulfate  fragments  in  the  AMS  spectrum  are 
m/z  48  (has  minor  organic  interference  based  on  unpublished 
high  resolution  data  of  a  similar  fuel  type),  m/z  64,  m/z  80, 
m/z  81  and  m/z  98.  The  ratios  of  m/z  64,  m/z  80,  m/z  81 
and  m/z  98  to  m/z  48  were  calculated  for  low  smoke  levels 
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Fig.  4.  Emission  ratio  plots  for  (a)  ACH3COOH/ACO  and  (b)  AC2  H4/ACO  measured  by  AFTIR.  (c)  AO  A/ A  CCA  measured  by  AMS  and 
NDIR.  (d)  ArBC/AC02  measured  by  SP2  and  NDIR. 


where  organic  interference  at  these  masses  is  low,  and  were 
then  applied  to  the  higher  loading  smoke  samples.  The  is¬ 
sues  related  to  calculating  inorganic  aerosol  concentrations 
from  UMR  AMS  data  during  episodes  of  high  organic  load¬ 
ing  is  an  active  area  of  research  and  will  be  presented  in  fu¬ 
ture  work. 

For  any  carbonaceous  fuel,  the  source  ERs  can  be  used  to 
calculate  emission  factors  (EF)  expressed  in  grams  of  com¬ 
pound  emitted  per  kilogram  of  biomass  burned  (on  a  dry 
weight  basis)  using  the  carbon  mass-balance  method  (as  ex¬ 
plained  in  Yokelson  et  al.,  1996,  1999).  The  carbon  mass 
balance  method  is  most  accurate  when  the  fraction  of  car¬ 
bon  in  the  fuel  is  precisely  known  and  all  the  burnt  carbon 
is  volatilized  and  detected.  Based  on  carbon  content  analysis 
of  about  50  chaparral  vegetation  samples  during  a  companion 
study  (Burling  et  al.,  2010)  and  literature  values  for  similar 
fuels  (Susott  et  al.,  1996;  Ebeling  and  Jenkins,  1985),  we 
assumed  a  carbon  content  of  50  %  by  mass  on  a  dry  weight 
basis  for  the  fuels  that  burned  in  the  Williams  Fire.  The  ac¬ 
tual  fuel  carbon  percentage  likely  varied  from  this  by  less 
than  a  few  percent.  EF  scale  linearly  with  the  assumed  fuel 
carbon  fraction.  Total  C  in  this  study  was  determined  from 
the  sum  of  the  C  from  AFTIR  species,  rBC,  and  an  organic 
carbon  (OC)  value  based  on  the  assumption  that  OA/OC  was 
1.6  in  fresh  smoke  (Yokelson  et  al.,  2009).  This  sum  could 
underestimate  the  total  carbon  by  1-2  %  due  to  unmeasured 
C.  That  would  lead  to  a  slight,  across  the  board  overestimate 
of  our  calculated  EFby  a  factor  of  1-2  %  (Akagi  et  al.,  201 1). 


3  Results  and  discussion 

3.1  Initial  emissions  compared  to  other  chaparral  fires 

The  fire-average  molar  emission  ratios  for  the  Williams  Fire 
(computed  as  described  in  Sect  2.4.2)  are  presented  in  Ta¬ 
ble  2.  A  detailed  comparison  of  the  Williams  Fire  initial 
emissions  (ER  and  EF)  to  the  emissions  from  other  chaparral 
fires  sampled  during  SLOBB  and  chaparral  fires  sampled  in 
other  studies  is  presented  by  Burling  et  al.  (2011).  Here  we 
simply  note  that  the  Williams  Fire  had  similar  ERs,  EFs,  and 
modified  combustion  efficiency  (MCE,  0.933)  to  those  re¬ 
ported  in  the  literature  for  other  chaparral  fires  (Radke  et  al., 
1991;  Hardy  et  al.,  1996;  Burling  et  al.,  2011).  The  one  ex¬ 
ception  is  that  Radke  et  al.  (1991)  reported  a  study-average 
EF(PM2.5)  (15.93  ±  6.88  g  kg-1)  that  is  significantly  higher 
than  our  EF(PM2.s)  or  the  literature  average  (Burling  et  al., 
2011,  ~8gkg_1),  which  likely  reflects  that  they  sampled 
more  of  the  smoldering  phase  near  the  end  of  the  fires,  as 
smoldering  normally  produces  higher  EF(PM)  (Ferek  et  al., 
1998).  Radke  et  al.  (1991)  note  that  their  spot  measurements 
of  particle  EF  increased  by  as  much  as  a  factor  of  three  as  the 
fraction  of  smoldering  combustion  increased.  Rather  than 
measure  samples  of  smoldering  combustion  as  the  fire  died 
down,  we  focused  first  on  characterizing  the  emissions  that 
were  produced  by  the  bulk  of  the  fuel  consumed  and  then  on 
the  evolution  of  those  emissions  in  the  downwind  plume. 
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3.2  Plume  evolution 

Many  chemical  and  physical  changes  can  affect  nascent 
smoke  as  it  begins  to  mix  with  ambient  air.  These  changes  in¬ 
clude  photochemical  transformations,  mixing  with  biogenic 
and/or  anthropogenic  pollutants,  cloud  processing,  coagu¬ 
lation,  and  gas-to-particle  conversion  (Reid  et  al.,  1998). 
The  complexity  of  the  aging  scenario  can  vary  greatly,  de¬ 
pending  on  factors  such  as  temperature,  dispersion  altitude, 
time  of  day,  humidity,  cloud  cover,  the  initial  emissions, 
and  the  proximity  of  other  emission  sources.  In  this  study, 
the  Williams  Fire  produced  an  isolated,  nearly  “undisturbed” 
plume,  which  provides  a  rare,  valuable  opportunity  to  test 
the  mechanism  of  a  photochemical  model.  A  companion 
paper  will  use  our  observations  of  the  Williams  Fire  to  test 
and  evaluate  an  advanced  smoke  plume  chemistry  model  (Al¬ 
varado  and  Prinn,  2009). 

The  evolution  of  the  reactive  species  we  could  measure 
over  the  course  of  4h  following  emission  is  summarized  in 
Table  2.  To  illustrate  photochemical  changes  in  the  plume 
it  is  useful  to  have  a  best  estimate  of  the  starting  NEMR 
and  uncertainty  to  compare  to  downwind  NEMR.  We  use  the 
slope -based  fire-average  ER  as  our  best  estimate  of  the  likely 
starting  NEMR  for  primary  species  measured  in  any  individ¬ 
ual  smoke  parcel  downwind.  Sloped  based  initial  NEMRs 
were  also  calculated  for  secondary  species  such  as  ozone 
and  PAN  (shown  in  Table  2)  though  the  concept  of  an  ER 
does  not  strictly  apply  to  these  species.  We  use  the  stan¬ 
dard  deviation  in  the  mean  of  the  ten  source  NEMRs  as  our 
best  estimate  of  the  uncertainty  in  the  starting  NEMR  for  all 
species  with  downwind  data  (also  shown  in  Table  2).  Stable 
species  that  did  not  change  significantly  for  4  h  after  emis¬ 
sion  and  species  that  dropped  below  our  detection  limit  in 
that  time  are  also  included  in  Table  2  just  to  provide  a  more 
complete  reference.  The  uncertainty  in  the  ER  shown  for 
these  species  is  the  standard  error  in  the  slope  and  no  down¬ 
wind  data  is  shown  except  for  rBC.  Absolute  uncertainties 
in  the  ERs  and  NEMRs  have  a  component  from  random  and 
systematic  instrumental  errors,  sampling  representativeness, 
and  other  issues;  but  many  of  these  errors  should  not  greatly 
impact  the  study  of  plume  evolution.  As  an  example  we  note 
that  downwind  smoke  was  only  sampled  in  the  high  altitude 
plume  and  the  low-altitude  source  samples  could  potentially 
represent  smoke  that  instead  detrained  in  the  low  elevation 
plume.  However,  the  source  ERs  were  not  significantly  de¬ 
pendent  on  altitude.  In  a  few  cases  we  will  illustrate  this  by 
comparing  the  downwind  NEMRs  to  an  alternate  ER  based 
solely  on  the  high  altitude  source  samples  and  show  that  the 
impact  of  this  was  small. 

At  this  point  it  should  be  clear  that  fires  do  not  produce 
“ideal”  plumes.  Instead  there  is  some  temporal  and  spatial 
fluctuation  about  the  average  ER  and  smoke  parcels  mix  with 
both  background  air  and  other  parts  of  the  plume,  which  ex¬ 
plains  some  of  the  variation  in  the  downwind  NEMRs.  We  fit 
the  evolution  of  the  NEMRs  to  various  “best  fit”  functions  to 
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Fig.  5.  AO3/ACO  vs.  time  since  emission  (h).  Individual  source 
samples  are  not  shown  and  we  instead  show  the  initial  negative 
AO3/ACO  NEMR  (solid  circle).  The  error  bar  reflecting  the  stan¬ 
dard  deviation  in  the  mean  of  the  ten  source  measurements  is  too 
small  to  be  visible.  The  downwind  error  bar  is  an  estimate  of  the 
uncertainty  in  the  oldest  individual  downwind  sample. 

estimate  the  final  NEMR  at  the  age  of  the  last  measurement 
(4-4.5  h  time  since  emission),  to  assist  in  discerning  differ¬ 
ences  between  the  final  downwind  NEMR  and  the  ER,  and 
to  estimate  lifetimes.  We  report  the  final  downwind  NEMR 
and  estimated  uncertainty  for  the  measured  reactive  species 
in  Table  2.  The  uncertainty  shown  for  the  final  NEMR  is  the 
uncertainty  in  individual  downwind  samples,  but  the  uncer¬ 
tainty  in  whether  a  trend  is  occurring  is  smaller  than  shown 
by  these  uncertainties.  We  note  that  the  functions  used  were 
selected  empirically  based  on  how  well  they  fit  the  data  and 
cannot  necessarily  be  rigorously  derived  from  the  complex, 
sometimes  unknown,  underlying  chemistry.  Nor  do  we  as¬ 
sume  that  these  fits  describe  the  behavior  beyond  the  time 
since  emission  measured  here.  We  show  relative  growth  fac¬ 
tors  that  compare  the  downwind  and  initial  NEMRs  (or  ERs) 
over  this  time  period  to  summarize  the  magnitude  of  species 
evolution  that  occurred  in  this  photochemically  active  plume. 
Detailed  discussions  of  these  changes  are  presented  in  the  re¬ 
mainder  of  this  section. 

3.2.1  Ozone 

Ozone  is  an  important  secondary  product  formed  in  biomass 
burning  smoke.  Figure  5  shows  the  AO3/ACO  NEMR 
tracked  over  the  course  of  estimated  times  since  emission  of 
Oh  to  ~4.5h.  The  initial  AO3/ACO  NEMR  (solid  circle) 
is  negative  because  the  background  air  contained  ~50ppb 
of  ozone  that  was  quickly  destroyed  by  fast  reaction  with 
NO  emitted  by  the  fire.  Negative  initial  AO3/ACO  NEMRs 
are  common  when  sampling  fresh  smoke  (Yokelson  et  al., 
2003a).  The  y-intercept  of  a  linear  trendline  was  forced  to  the 
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Table  2.  Emission  ratios  (ER),  initial  and  downwind  NEMR  for  the  Williams  Fire. 


ER  (mol  mol  *) 
or  initial  NEMR 

Uncertainty11 

Downwind  NEMR 
(4^1.5  h  since  emission) 

Uncertainty*5 

Factor  Change 
after  4^1.5  h 

Uncertainty 

Ratioed  to  ACCU 

CO 

7.13  x  10“2 

5.50  x  10“3 

Ratioed  to  ACO 

Methane  (CH4) 

8.72  x  10"2 

2.40  x  10“3 

Acetylene  (C2H2) 

2.72  x  lO"3 

2.10  x  10“4 

Propylene  (C3H6) 

4.80  x  10“ 3 

8.00  x  10”5 

Phenol  (C6H5OH) 

1.48  x  10“ 3 

1.70  x  10“4 

Furan  (C4H4O) 

1.46  x  10“3 

8.00  x  10“5 

Methanol  (CH3OH) 

1.67  x  10“2 

7.00  x  10“4 

Hydrogen  Cyanide  (HCN) 

1.30  x  10“2 

1.10  x  10“3 

Nitrous  Acid  (HONO) 

4.00  x  10“ 3 

6.10  x  10“4 

Ethylene  (C2H4) 

1.28  x  10“2 

1.10  x  10”3 

6.01  x  10“3 

7.15  x  10“4 

0.47 

0.07 

Glycolaldehyde  (HCOCH2OH) 

1.87  x  10“5 

1.18  x  10“5 

2.52  x  10“3 

2.12  x  10“ 3 

135 

142 

Formic  Acid  (HCOOH) 

6.59  x  10“4 

1.53  x  10“4 

4.84  x  10“3 

1.65  x  10“3 

7.34 

3.03 

Acetic  Acid  (CH3COOH) 

1.41  x  10“2 

1.19  x  10”3 

2.45  x  10“2 

5.69  x  10“3 

1.73 

0.43 

Formaldehyde  (HCHO) 

1.65  x  10“2 

2.59  x  10“3 

2.42  x  10“2 

4.74  x  10“3 

1.47 

0.37 

Ammonia  (NH3 ) 

3.84  x  10“2 

9.69  x  10“3 

1.91  x  10“2 

2.65  x  10“3 

0.50 

0.14 

Ozone  (O3) 

-5.13  x  10“3 

1.13  x  10“3 

1.02  x  10“* 

2.16  x  10“2 

N/A 

N/A 

Ratioed  to  ACO2 

Peroxy  Acetyl  Nitrate  (PAN) 

4.75  x  10“5 

2.04  x  10“5 

5.10  x  10“4 

1.21  x  10“4 

10.74 

5.27 

Nitric  Oxide  (NO) 

8.21  x  10“4 

1.17  x  10“4 

Nitrogen  Dioxide  (NO2) 

1.64  x  10“3 

9.00  x  10“5 

NOx  (as  NO) 

2.30  x  10“ 3 

3.62  x  10“4 

4.60  x  10“4 

2.30  x  10“4 

0.20 

0.10 

PM2.5  (light  scattering)0,11 

5.15  x  10“ 3 

1.33  x  10“3 

1.29  x  10“2 

1.89  x  10“3 

2.50 

0.74 

OCc 

2.22  x  10“3 

5.36  x  10“4 

OA° 

3.55  x  10“ 3 

8.57  x  10“4 

2.83  x  10“3 

2.52  x  10“4 

0.80 

0.21 

Sulfate  (SO2-)0 

3.58  x  10“6 

3.28  x  10“6 

6.72  x  10“6 

2.42  x  10“6 

1.88 

1.85 

Nitrate  (NOj")c 

9.61  x  10“5 

7.09  x  10“5 

3.87  x  10“4 

1.39  x  10“4 

4.03 

3.30 

Ammonium  (NH|)° 

6.00  x  10“5 

3.95  x  10“5 

2.34  x  10“4 

8.42  x  10“5 

3.90 

2.93 

Chloride  (Cl-)0 

4.97  x  10“5 

2.15  x  10“5 

1.33  x  10“5 

4.79  x  10“6 

0.27 

0.15 

rBC0 

7.83  x  10“4 

3.48  x  10“4 

7.30  x  10“4 

1.51  x  10“4 

0.93 

0.46 

fTC  rBCe 

0.54 

0.07 

0.85 

0.01 

1.57 

0.20 

a  Error  calculated  as  the  standard  error  in  the  slope  when  no  downwind  NEMR  data  shown.  Otherwise,  uncertainty  was  calculated  as  the  standard  deviation  in  the  mean  of  the  source 
measurements  collected. 

k  Downwind  errors  shown  are  an  estimate  of  the  uncertainty  in  individual  downwind  samples.  The  uncertainty  in  whether  a  trend  is  occurring  is  smaller. 
c  Reported  as  a  mass  ratio. 

^  Converting  light  scattering  to  PM2.5  with  constant  mass  scattering  efficiency  not  valid  for  aged  smoke,  see  Sect.  3.3. 
e  Fraction  of  “thickly  coated"  rBC  particles.  All  numbers  reported  as  fractions. 


initial  slope-based  AO3/ACO  NEMR  at  time  t  —  0  h  and  the 
individual  source  measurements  are  not  shown.  The  down¬ 
wind  error  bar  is  an  estimate  of  the  uncertainty  in  the  old¬ 
est  individual  downwind  sample.  As  shown,  AO3/ACO  in¬ 
creased  from  —5.13  (±1.13)  x  10“  3  to  10.2  (±2.16)  x  10“2 
in  ~4.5h  since  smoke  emission.  Ozone  production  in  a 
temperate  forest  BB  plume  in  the  Pacific  Northwest  of  the 
US  was  noted  previously  by  Hobbs  et  al.  (1996),  who  re¬ 
ported  that  AO3/ACO  reached  1.5  x  10“2  in  30  min.  Akagi 
et  al.  (2011)  assessed  the  literature  on  O3  measurements  in 
BB  plumes  and  concluded  that  O3  formation  is  ubiquitous  in 


tropical  BB  plumes,  but  O3  formation  or  destruction  can  oc¬ 
cur  in  extratropical  BB  plumes.  Global  models  consistently 
conclude  that  BB  is  a  major  global  source  of  tropospheric  O3 
(Fishman  et  al.,  2003;  Sudo  and  Akimoto,  2007).  However, 
detailed  plume  scale  studies  find  that  reproducing  detailed 
observations  of  O3  formation  is  difficult  (Trentmann  et  al., 
2005;  Alvarado  and  Prinn,  2009).  The  data  we  report  here 
provides  a  rare  opportunity  to  test  the  photochemical  mech¬ 
anism  in  a  model  with  low  uncertainty  in  transport  (M.  Al¬ 
varado,  personal  communication,  2012). 
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Fig.  6.  Ln(Aethylene/ACO)  vs.  time  since  emission  (h)  x  (kethylene 
—kQo)-  Average  OH  in  the  plume  from  ~10:30-15:16LT  is 
initially  estimated  from  the  absolute  value  of  the  slope  (5.70  ± 
1.05  x  106  molec  cm-3)  and  then  corrected  for  O3  oxidation  (see 
Sect.  3.2.2). 

3.2.2  Alkenes  and  estimated  OH 

Propylene  mixing  ratios  dropped  to  levels  near  our  detection 
limits  within  ~30  min  due  to  oxidation  and  dilution.  The 
low  signal-to-noise  for  C3H6  in  the  aged  smoke  made  it  dif¬ 
ficult  to  quantify  the  decay  rate.  We  also  observed  a  rapid, 
but  readily  quantifiable,  decrease  in  AC2H4/ACO  (Fig.  6). 
We  estimate  the  average  OH  concentration  in  the  plume  by 
considering  both  main  oxidation  processes  for  C2H4.  Re¬ 
action  with  OH  is  the  primary  removal  process  (kethyiene  = 
8.32x  10-12  cm3  molec-1  s-1;  Sander  et  al.,  2006).  Ozone 
oxidation  is  a  competing  reaction,  but  the  loss  of  C2H4 
is  relatively  small  (~7.6%)  given  the  lower  rate  constant 
(^ozone  =  1.7  x  10-18  cm3  molec-1  s-1;  Sander  et  al.,  2006) 
and  the  typical  O3  levels  (~100ppb)  measured  in  the  plume. 
As  a  first  step,  following  the  procedure  derived  by  Hobbs  et 
al.  (2003),  we  fit  a  line  to  the  decay  of  ln(AC2H4/ACO)  ver¬ 
sus  (time  since  emission)  x  (kethylene  —  kcoX  where  kco  = 
1.74  x  10-13  cm3  molec-1  s-1  (Sander  et  al.,  2006)  (Fig.  6). 
The  slope  in  the  plot  implies  an  average  OH  concentration 
of  5.70  (±1.05)  x  106  molec  cm-3.  We  reduce  this  average 
OH  by  7.6%  to  account  for  the  ethylene  lost  via  the  O3 
oxidation  channel  resulting  in  a  corrected  OH  estimate  of 
5.27  (±0.97)  x  106  molec  cm-3.  If  we  use  only  the  high- 
altitude  source  samples  it  raises  the  initial  AC2H4/ACO  by 
~12  %  and  that  would  also  raise  the  OH  estimate.  The  OH 
and  O3  were  unlikely  to  be  constant  in  the  plume  and  the 
background  CO  was  as  high  as  22  %  of  total  CO  in  the  most 
dilute  sample.  However,  omitting  the  2-4  most  dilute  sam¬ 
ples  or  the  oldest  samples  had  only  a  small  effect  on  the  OH 
retrieval.  Thus,  we  use  the  O3  corrected  OH  concentration  of 
5.27  (±0.97)  x  106  molec  cm-3  as  a  reasonable  estimate  of 
average  OH  in  this  plume. 
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Fig.  7.  Increase  in  ACH3COOH/ACO  (blue)  and  AHCOOH/ACO 
(green)  vs.  time  since  emission  (h).  Individual  source  samples  are 
not  shown  and  we  instead  show  ERs  (solid  blue  and  solid  green 
circles,  respectively)  with  an  error  bar  reflecting  the  standard  devia¬ 
tion  in  the  mean  of  the  ten  source  measurements.  The  ER(HCOOH) 
error  bar  is  too  small  to  be  visible.  Downwind  error  bars  are  an  esti¬ 
mate  of  the  uncertainty  in  the  oldest  individual  downwind  samples. 

The  OH  levels  largely  govern  the  lifetime  of  VOCs  and 
semi-volatile  organic  compounds  (SVOCs)  and  the  rate  of 
formation  of  soluble  or  less  volatile  secondary  products  (Se¬ 
infeld  and  Pandis,  2006).  Other  oxidants  such  as  HO2  and 
RO2  are  also  likely  elevated  in  high  OH  environments  (Ma¬ 
son  et  al.,  2001).  Similarly  high  concentrations  of  OH,  some 
greater  than  107  molec  cm-3,  have  been  inferred  or  measured 
in  African  and  Mexican  biomass  burning  plumes  within  the 
first  hour  of  plume  processing  (Hobbs  et  al.,  2003;  Yokel- 
son  et  al.,  2009).  The  Williams  Fire  OH  was  likely  lower 
than  that  in  the  African  and  Mexican  plumes  largely  be¬ 
cause  of  less  intense  solar  radiation  that  would  occur  in  mid- 
November  at  ~35°  N  (Tables  3.7  and  3.9  in  Finlayson-Pitts 
and  Pitts,  2000).  However,  the  Williams  Fire  OH  is  by 
no  means  depleted  and  likely  contributes  to  the  O3  produc¬ 
tion  described  above  and  many  of  the  additional  rapid  post¬ 
emission  changes  described  in  the  following  sections. 

3.2.3  Organic  acids  and  aldehydes 

The  aging  plots  for  AFTIR  and  AMS  species  from  this  point 
forward  (Fig.  7-11)  show  the  ER  of  A  Y/ACO  or  A  A/ACO2 
at  an  estimated  time  since  emission  of  t  =  0  h  with  the  y- 
intercept  of  a  best-fit  trendline  forced  to  the  initial  slope- 
based,  fire-average  ER,  shown  in  Table  2.  The  uncertainty 
in  the  ER  is  an  estimate  of  source  variability  from  the  1  a 
standard  deviation  in  AY/ AGO  in  the  ten  source  samples 
(described  in  Sect.  3.2).  The  uncertainty  shown  for  the  down¬ 
wind  samples  is  an  estimate  of  the  uncertainty  in  individual 
downwind  samples. 
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We  measured  formation  of  both  acetic  acid  and  formic  acid 
over  the  course  of  4.5  h  (Fig.  7).  The  ACH3COOH/ACO  ra¬ 
tio  increased  from  1.41  (±0.119)  x  10-2  to  2.45  (±0.569) 
x  10“2  over  4.5  h,  a  factor  of  -1.73  ±0.43  (Table  2). 
Over  the  same  aging  period  AHCOOH/ACO  increased  from 
6.59  (±1.53)  x  10“4  to  4.84  (±1.65)  x  10“3,  a  factor  of 
7.34  ±3.03.  Organic  acids  are  a  major  component  of 
biomass  burning  plumes  and  their  post-emission  behavior  in 
BB  plumes  is  summarized  in  Yokelson  et  al.  (2009).  The  in¬ 
crease  in  acetic  acid  in  the  Williams  Fire  plume  was  slower 
than  the  rapid  increase  in  ACH3COOH/ACO  to  —0.09  in 
less  than  1  h  seen  in  African  savanna  fire  plumes  (Yokelson  et 
al.,  2003a).  However,  formic  acid  did  not  increase  downwind 
of  those  savanna  fires.  In  the  Williams  Fire  plume,  and  many 
other  studies,  the  post-emission  increases  in  organic  acids 
were  greater  than  the  amount  of  known  precursors,  which 
suggests  that  organic  acids  may  be  the  oxidation  products 
of  currently  unidentified  co-emitted  NMOC  (Yokelson  et  al., 
2009).  At  higher  RH,  heterogeneous  reactions  may  also  be  a 
source  of  formic  acid  (Tabazadeh  et  al.,  1998). 

The  glycolaldehyde  NEMR  to  ACO  increases  by  a  small 
total  amount  downwind,  but  with  a  high  uncertainty  (Ta¬ 
ble  2).  Glycolaldehyde  is  a  major  product  of  cellulose  py¬ 
rolysis  (Richards,  1987;  Radlein  et  al.,  1991)  and  smolder¬ 
ing  combustion  (Yokelson  et  al.,  1997),  and  also  a  known 
byproduct  of  hydroxyl  initiated  oxidation  of  ethylene  and 
isoprene.  Being  highly  water  soluble,  aqueous-phase  reac¬ 
tions  of  glycolaldehyde  may  be  important  in  the  formation 
of  secondary  organic  aerosol  (Perri  et  al.,  2009)  via  the  for¬ 
mation  of  lower  volatility  compounds.  Formaldehyde  can 
be  produced  post-emission  in  plumes  from  the  oxidation  of 
NMOC,  but  also  lost  by  photolysis  and  reaction  with  OH  or 
HO2.  HCHO  photolysis  is  an  important  source  of  OH  (Ma¬ 
son  et  al.,  2001).  In  the  Williams  Fire  plume,  AHCHO/ACO 
increased  during  the  first  hour  after  emission  and  then  leveled 
off  (Table  2,  graph  not  shown). 

3.2.4  Gas-particle  conversion  and  formation  of 

reservoir  species:  nitrogen-containing  compounds 

Fitting  an  exponential  decay  to  the  molar  ratio  of 
ANH3/ACO  (not  shown)  implies  that  roughly  half  of  the  ini¬ 
tial  NH3  was  lost  in  —4.5  h.  Much  of  the  decrease  between 
the  downwind  NEMRs  and  the  ER  occurs  in  the  first  hour 
following  emission,  suggesting  rapid,  initial  NH3  removal 
by  chemical  and  physical  aging.  In  Fig.  8  we  compare  the 
decrease  in  the  molar  ratio  of  ANH3/ACO2  (derived  from 
AFTIR  data)  to  the  molar  increase  in  ANHJ/ACO2  (mea¬ 
sured  by  the  AMS  and  NDIR)  over  4h,  which  is  the  ap¬ 
proximate  time  since  emission  of  the  oldest  AMS  sample. 
The  fit  implies  a  significant  decay  in  ANH3/ACO2  along 
with  a  clear  increase  in  ANH^/ACCb  (Fig.  8a  and  b,  re¬ 
spectively).  We  note  the  apparent  high  variability  in  the 
source  measurements  for  NH3  and  emphasize  that  this  is 
an  artifact  arising  from  compression  to  “time  zero”  of  all 
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Fig.  8.  (a)  ANH3/ACO2  molar  ratio  vs.  time  since  emission  (h). 
The  decay  in  ANH3/ACO  vs.  time  since  emission  is  slightly  better 
correlated,  but  we  show  ANH3/ACO2  here  for  direct  comparison 
to  the  AMS  data  in  Fig.  8b.  Points  were  fit  to  an  exponential  trend- 
line.  (b)  ANH^"/AC02  molar  ratio  vs.  time  since  emission  (h). 
Note  difference  in  y-scale.  Points  were  fit  to  a  logarithmic  trend- 
line.  For  both  (a)  and  (b),  ERs  (solid  circles)  are  shown  with  an 
error  bar  reflecting  the  standard  deviation  in  the  mean  of  the  ten 
source  measurements.  We  retain  all  source  samples  to  show  the  in¬ 
herent  variability  in  source  measurements.  Downwind  error  bars 
are  an  estimate  of  the  uncertainty  in  the  oldest  individual  downwind 
samples. 

the  ANH3/ACO2  ER  samples,  which  increased  in  mono¬ 
tonic  fashion  from  the  beginning  to  end  of  the  source  sam¬ 
pling  period.  The  ANHJ/ACO2  ratio  increases  post  emis¬ 
sion  by  a  factor  of  3.90  ±2.93  (Table  2),  with  this  growth 
accounting  for  a  fraction  of  0.36  ±0.31  of  the  gaseous  am¬ 
monia  lost  on  a  molar  basis.  NH3  reacts  quickly  with  HNO3 
to  form  solid  NH4NO3  and  NH3  may  also  react  with  or¬ 
ganic  acids  to  form  other  ammonium  compounds  in  both  the 
solid  and  liquid  phases  (though  the  low  RH  in  the  Williams 
Fire  plume  favors  reactions  in  the  solid  phase,  see  below). 
However,  since  only  about  one-third  of  the  ammonia  lost 
was  detected  as  ammonium,  other  loss  processes  must  be 
considered.  Ammonia  removal  pathways  competing  with 
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ammonium  formation  include  the  oxidation  of  NH3  by  OH 
(Dentener  and  Crutzen,  1994),  but  this  reaction  proceeds  at 
a  relatively  slow  rate  at  typical  plume  concentrations  ( k  = 
1.6  x  10-13  cm3  molec-1  s,  Sander  et  al.,  2006)  and  would 
not  significantly  reduce  the  amount  of  free  ammonia  over 
the  course  of  4  h.  No  precipitation  or  even  clouds  impacted 
the  plume  so  wet  deposition  was  not  a  factor,  but  a  plausi¬ 
ble  loss  process  was  the  dry  deposition  of  NH3  since  it  is  a 
“sticky”  molecule  known  to  adhere  on  surface  sites  (Yokel- 
son  et  al.,  2003b)  and  undergo  dry  deposition  readily  (Sein¬ 
feld  and  Pandis,  2006).  Dry  deposition  of  NH3  to  any  water 
coating  that  could  be  present  on  the  surface  of  aerosol  par¬ 
ticles  can  remove  ammonia  from  the  gaseous  state  (Asman 
et  al.,  1998),  however,  over  99%  of  gaseous  ammonia  ab¬ 
sorbing  into  the  liquid  phase  would  exist  as  ammonium  and 
be  detected  by  the  AMS.  Dry  particles  can  also  provide  sur¬ 
face  sites  for  ammonia  adsorption  to  occur.  However,  NH3 
adsorbed  to  particles  would  also  likely  be  detected  by  the 
AMS  as  ammonium,  since  particulate  ammonium  fragments 
to  m/z  17,  16,  and  15,  following  vaporization  and  ionization. 
Thus,  it  is  not  likely  that  dry  deposition  accounts  for  the  ma¬ 
jority  of  the  observed  NH3  loss. 

Other  plausible  fates  of  NH3  may  have  involved  the  for¬ 
mation  of  new  ammonium  nitrate  particles  at  the  plume  tem¬ 
peratures  of  12-15  °C  or  condensation  of  ammonium  nitrate 
onto  existing  particles  that,  in  either  case,  are  outside  the 
most  efficient  detection  range  of  the  AMS  (50-800  nm).  Due 
to  the  low  RH  environment  of  the  Williams  Fire,  nascent 
particles  were  probably  not  water  coated  and  grew  very 
slowly.  While  Semuniuk  et  al.  (2007)  observed  rapid  hy¬ 
groscopic  growth  of  fresh  smoke  particles  in  the  range  of 
55-100  %RH,  the  Williams  Fire  RH  range  of  -11-26%  is 
much  lower,  and  also  lower  than  the  deliquescence  and  efflo¬ 
rescence  points  of  most  common  salts  (>60%)  (Tang  and 
Munkelwitz,  1993,  1994a,  b).  Thus,  new  particle  forma¬ 
tion  and  condensation  on  those  particles  could  have  been 
occurring  without  being  detected  if  particles  were  growing 
slow  enough  to  evade  AMS  detection  even  after  several  hours 
since  emission.  We  also  can  not  discount  the  possibility  that 
ammonium  nitrate  was  incorporated  in  particles  that  were 
larger  than  800  nm. 

We  roughly  compare  the  Williams  Fire  observations  with 
those  in  a  Yucatan  BB  plume  in  a  70-100  %  RH  environment 
where  ammonium  formation  accounted  for  roughly  32  %  of 
the  estimated  starting  NH3  in  just  1.4  h  (Yokelson  et  al., 
2009).  Ammonium  formation  accounted  for  only  —13% 
of  the  measured  initial  NH3  over  1 .4  h  in  the  Williams  Fire 
plume.  Faster  ammonium  formation  accounting  for  more 
of  the  initial  ammonia  may  have  occurred  in  the  Yucatan 
due  to  increased  deposition  to  wetter  particles,  faster  particle 
growth,  or  faster  HNO3  formation  under  tropical  photochem¬ 
ical  conditions. 

Figure  9  shows  the  decrease  in  ANOx/AC02  (blue)  and 
corresponding  increases  in  both  ANOjVACCH  (green)  and 
APAN/ACO2  (light  blue)  during  4h  after  emission.  Pro¬ 
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Fig.  9.  ANOx  (blue),  APAN  (light  blue),  and  ANO^"  (green)  ra- 
tioed  to  ACO2  (molar)  vs.  time  since  emission  (h).  Corresponding 
emission  ratios  (ER)  are  shown  as  solid  blue,  light  blue,  and  green 
circles,  respectively.  ER  error  bars  show  the  standard  deviation  in 
the  mean  of  the  ten  source  measurements.  The  ER(PAN)  error  bar 
is  too  small  to  be  visible.  Downwind  error  bars  are  an  estimate  of 
the  uncertainty  in  oldest  individual  downwind  samples.  Points  were 
fit  to  logarithmic  trendlines. 


cesses  that  remove  NOx  from  the  plume  or  transfer  it  to 
reservoir  species  such  as  PAN  or  alkyl  nitrates  are  impor¬ 
tant,  as  NOx  is  a  precursor  for  O3  formation.  A  major  NOx 
loss  process  is  the  association  reaction  of  OH  with  NO2  to 
form  HNO3,  which  can  condense  to  form  particle  nitrate.  At 
the  average  temperature  and  pressure  of  the  Williams  Fire 
plume,  the  NO2  association  reaction  occurs  with  a  pseudo 
second  order  rate  constant  of  1.04  x  10“ 11  cm3  molec-1  s-1 
(Sander  et  al.,  2006)  and  a  pseudo  first  order  rate  con¬ 
stant  (using  our  derived  average  OH  concentration  from 
Sect.  3.2.2)  of  5.47  x  10-5  s-1  for  an  NO2  lifetime  of 
—5.1  h.  Thus,  in  the  4  h  that  we  monitored  the  plume  this  re¬ 
action  would  be  expected  to  convert  roughly  54  %  of  the  ini¬ 
tial  NOx  to  HNO3  (and  potentially  nitrate).  Gas-phase  nitric 
acid  was  not  measured  in  the  Williams  Fire  plume,  however 
other  studies  have  reported  that  there  was  no  measureable 
enhancement  of  gas-phase  HNO3  in  any  of  the  BB  plumes 
measured  (Yokelson  et  al.,  2009;  Alvarado  et  al.,  2010).  The 
reason  for  this  is  likely  because  HNO3  can  react  rapidly  with 
NH3  to  form  NH4NO3  (as  discussed  above),  and  HNO3  (or 
NH4NO3)  have  high  deposition  velocities  to  wet  and  dry  sur¬ 
faces.  In  the  Williams  Fire,  all  these  processes  could  theoret¬ 
ically  remove  HNO3  from  the  plume,  but  loss  by  deposition 
to  wet  surfaces  is  unlikely  due  to  the  low  RH. 

The  ANO^/ACCH  molar  ratio  was  initially  very  low,  but 
secondary  formation  of  nitrate  occurred  rapidly  as  can  be 
seen  from  the  logarithmic  trendline  in  Fig.  9  (green).  The 
trend  in  ANCfjVACCH  is  similar  to  that  of  ANH^/ACCH, 
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suggesting  an  initial  burst  in  these  inorganic  particle  species 
in  the  first  hour  since  emission.  After  4h,  AMS-detected 
nitrate  accounted  for  only  ~22  %  of  the  HNO3  (molar  ba¬ 
sis)  predicted  to  form  based  on  our  derived  average  OH  con¬ 
centration.  The  unaccounted  for  HNO3  comprises  ~42  % 
of  the  initial  NOx.  Some  of  the  missing  HNO3  could  have 
been  present  in  the  plume  as  nitrate  in  particles  that  disso¬ 
ciated  in  the  AMS  inlet.  However,  inlet  losses  in  the  AMS 
due  to  evaporation  probably  did  not  significantly  affect  our 
measurements  given  the  modest  temperature  increase  within 
the  aircraft  (3-1 1  °C)  and  the  inlet  residence  time  of  ~2.5  s. 
Huffman  et  al.  (2009)  outfitted  their  airborne  AMS  inlet  with 
a  thermodenuder  and  measured  nitrate  evaporation  rates  of 
~0.2%C-1s-1.  As  mentioned  above,  it  is  plausible  that 
NH3  and  HNO3  formed  new  (or  condensed  on  existing)  par¬ 
ticles  that  were  outside  the  AMS  size  detection  range.  Miss¬ 
ing  HNO3  ratioed  to  missing  ammonium  (~1.3)  does  sug¬ 
gest  that  small  or  large  particles  containing  ammonium  ni¬ 
trate  could  account  for  substantial  undetected  signal  due  to 
both  ions,  but  that  some  gas  phase  HNO3  may  have  also  been 
present  that  we  were  unable  to  measure  with  the  instruments 
onboard  the  Twin  Otter.  The  possible  importance  of  ammo¬ 
nium  nitrate  in  small  particles  in  fresh  BB  plumes  is  further 
supported  by  the  observations  of  Yokelson  et  al.  (2009).  In 
that  study,  49  ±  16%  of  HNO3  expected  to  form  from  the 
reaction  of  NO2  and  OH  was  observed  as  nitrate  in  a  high 
RH  BB  plume  after  1.4  h  since  emission.  Further,  HNO3 
was  measured,  but  not  elevated  above  background  in  the 
plume.  Thus,  in  that  high  RH  environment  (70-100%),  par¬ 
ticle  growth  could  have  been  approximately  six  times  faster 
than  in  the  Williams  Fire,  but  still  slow  enough  that  half 
the  nitrate  mass  was  in  particles  too  small  to  detect  with  the 
AMS. 

Formation  of  reservoir  species  such  as  PAN  also  causes 
NOx  to  decrease  from  its  initial  level  (Seinfeld  and  Pandis, 
2006;  Alvarado  et  al.,  2010).  PAN  ratioed  to  CO2  shows 
rapid  initial  growth  (Fig.  9,  light  blue)  similar  to  nitrate  and 
ammonium.  Overall,  APAN/ACO2  increases  by  a  factor  of 
10.7  ±5.3  in  4h  since  emission  (Table  2)  and  accounts  for 
22  ±  6  %  of  initial  NOx  on  a  molar  basis  after  4  h.  Thus,  the 
total  NOx  decrease  due  to  the  observed  PAN  formation  plus 
the  calculated  NOx  loss  via  reaction  with  OH  is  76  %,  which 
agrees  remarkably  well  with  the  observed  total  NOx  decrease 
of  '-80  %  seen  in  Fig.  9.  However,  the  uncertainties  are  high 
enough  that  formation  of  other  reservoir  species  such  as  gas- 
phase  HNO3,  alkyl  nitrates,  or  other  peroxy  nitrates  could 
also  be  significant.  The  decrease  in  ANOx/AC02  over  4h 
shows  how  NOt,  which  is  available  for  photolysis  and  O3 
formation  in  the  fresh  plume,  was  quickly  converted  to  PAN, 
which  can  then  be  transported  considerable  distances  before 
thermally  decomposing  and  releasing  NO2  for  O3  produc¬ 
tion  far  from  the  source  (Val  Martin  et  al.,  2006).  After  4h 
since  emission,  the  excess  PAN  is  roughly  double  the  excess 
NOT.  Alvarado  et  al.  (2010)  also  observed  a  PAN/nitrate 
ratio  of  2  in  boreal  forest  fire  plumes  of  qualitatively  sim- 
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ilar  time  since  emission.  In  the  case  of  the  Williams  Fire, 
we  are  hesitant  to  interpret  the  PAN/nitrate  ratio  as  a  mea¬ 
sure  of  the  APAN/AHNO3  branching  ratio  since  the  nitrate 
fails  to  account  for  all  the  NOx  calculated  to  be  converted 
to  HNO3.  Smoke  plume  simulation  studies  by  Mason  et 
al.  (2001)  found  that  PAN  formation  becomes  increasingly 
dominant  over  HNO3  formation  further  downwind.  While 
we  cannot  account  for  all  NOx  lost  to  the  HNO3  channel,  we 
do  observe  that  the  APAN/ANO^  ratio  grew  by  a  factor  of 
~2  over  the  course  of  4  h  since  emission. 

3.2.5  Evolution  of  sulfate  and  chloride 

The  AS04  /ACO2  mass  ratio  steadily  increased  nearly  two¬ 
fold  from  3.58  (±3.28)  x  10“6  to  6.72  (±2.42)  x  10“6  in 
4h  since  emission  (Table  2,  Fig.  10a).  The  high  uncer¬ 
tainty  reported  is  due  to  the  large  variability  in  the  source 
measurements,  but  a  clear  increasing  trend  with  time  since 
emission  is  evident.  Unlike  the  rapid  initial  formation  of 
nitrate  and  ammonium,  the  sulfate  growth  is  gradual.  This 
is  expected  since  OH  reacts  relatively  slowly  with  SO2 
(8.8  x  10“ 13  cm3  molec-1  s_1,  Sander  et  al.,  2006),  and  the 
plume  was  too  dry  for  aqueous  sulfate  formation  to  be  sig¬ 
nificant.  In  a  Yucatan  BB  plume  characterized  by  high  initial 
emissions  of  SO2  and  high  RH,  Yokelson  et  al.  (2009)  ob¬ 
served  that  AS04_/ACO  increased  by  a  factor  of  8.8  ±2.7 
in  just  1 .4  h.  This  suggests  that  the  higher  RH,  OH,  and  ini¬ 
tial  SO2  contributed  to  much  faster  sulfate  formation  in  the 
Yucatan  plume.  In  addition,  Yokelson  et  al.  (2009)  observed 
ANHj"  to  be  about  equal  to  the  sum  of  twice  the  AS04  plus 
the  ANOij  (on  a  molar  basis)  in  that  plume.  Additionally, 
sulfate  had  accounted  for  ~60  %  as  much  negative  charge  as 
nitrate.  For  the  Williams  Fire,  a  graph  of  A  NOT  vs.  ANH4 
on  a  molar  basis  (Fig.  10b)  yields  a  slope  of  0.466  ±0.010 
(R2  —  0.95),  which  is  almost  identical  to  the  ratio  observed  in 
the  Yucatan  plume,  but  the  Williams  Fire  ANOjT  to  AS04_ 
molar  ratio  was  >100.  This  suggests  sulfate  was  not  a  sig¬ 
nificant  counter-ion  for  ammonium  in  this  plume  and  that  the 
particle  ammonium  was  pairing  with  negative  ions  other  than 
nitrate  and  sulfate. 

The  ACU/ACO2  ratio  decreases  with  time  since  emission 
in  a  well  correlated  fashion  (Fig.  10c),  except  for  one  sample 
collected  during  Flight  1 .  Chloride  may  decrease  because  it 
is  replaced  by  other  anions  such  as  nitrate  as  they  form  (Pos- 
fai  et  al.,  2003;  Li  et  al.,  2003)  although  nitrate  formation 
greatly  exceeds  chloride  losses.  The  fate  of  the  chloride  re¬ 
leased  from  the  particles  is  unknown,  but  it  is  unlikely  that 
most  of  the  released  chloride  constitutes  a  rapid  source  of 
Cl  atoms.  The  Cl  atom  has  an  appreciable  rate  coefficient 
with  C2H2  (Tabazadeh  et  al.,  2004),  and  we  did  not  observe 
a  change  in  AC2H2/ACO  over  4.5  h  since  emission.  In  any 
case,  the  amount  of  AC1-  is  insufficient  to  pair  with  a  signif¬ 
icant  fraction  of  the  ANHjJ"  during  the  later  stages  of  plume 
evolution  (Fig.  lOd)  and  the  total  charge  from  measured  inor¬ 
ganic  anions  (Cl“  +  twice  sulfate  +  nitrate)  is  only  sufficient 
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Fig.  10.  (a)  ASO^  /ACO2  mass  ratio  vs.  time  since  emission  (h).  (b)  ANO^  /ACO2  vs.  ANH^7AC02  on  a  molar  basis,  (c)  AC1  /ACO2 
mass  ratio  vs.  time  since  emission  (h)  fit  to  a  power  curve,  (d)  AC1-/ACC>2  vs.  ANH^ /ACO2  on  a  molar  basis.  In  both  (a)  and  (c)  individual 
source  samples  are  not  shown  and  we  instead  show  ERs  (solid  circles)  with  an  error  bar  that  is  the  standard  deviation  in  the  mean  of  the  ten 
source  measurements.  Downwind  error  bars  are  an  estimate  of  the  uncertainty  in  oldest  individual  downwind  samples. 


to  pair  with  ~50  %  of  the  ammonium  suggesting  that  the  lat¬ 
ter  may  also  be  pairing  with  unmeasured  organic  anions  such 
as  formate  or  acetate.  The  possible  significance  of  organic 
anions  is  further  supported  by  the  fact  that  ammonium  is  not 
the  only  significant  cation  in  BB  plumes.  For  instance  high 
levels  of  K+  are  normally  reported  in  BB  plumes  and  it  is 
often  used  as  a  tracer  for  BB  (Zhang  et  al.,  2010). 

3.3  Evolution  of  light  scattering  and  organic  aerosol 

The  trends  in  light  scattering  and  organic  aerosol  (OA)  were 
similar  when  we  used  source  ER  based  on  all  the  source 
samples  or  only  those  obtained  at  high  altitudes.  In  general, 
A/?scat  lor  Alight  scatteringVACCH  increased  by  a  factor  of 
2.50  ±0.74  (Table  2)  and  AZ?scat/ACO  increased  by  a  sim¬ 
ilar  factor  of  2.45  ±  0.55  over  the  course  of  4  h.  If  we  as¬ 
sume  a  constant  conversion  factor  with  time  since  emission 
(208  800  pg  sm~2,  Sect.  2.1.4)  between  scattering  (at  530  nm 
wavelength)  and  PM2.5,  then  this  implies  a  factor  of  ~2.50 
increase  in  APM2.5  (ratioed  to  ACO2)  as  shown  in  Fig.  11. 
This  increase  is  comparable  to  the  2.6  ±  1.5  growth  factor 
in  Alight  scattering/ ACO  seen  in  ~1.4h  in  an  isolated  BB 
plume  in  the  Yucatan  (Yokelson  et  al.,  2009).  In  the  case 
of  the  Yucatan  plume,  the  particle  growth  inferred  from  in¬ 
creased  scattering  (derived  using  the  same  methodology  as 


above)  was  confirmed  by  a  growth  factor  of  ~2.4  in  the  sum 
of  species  measured  by  AMS  ratioed  to  ACO.  The  growth  in 
AMS  species  was  driven  mostly  by  significant  secondary  OA 
(SOA)  formation  in  the  Yucatan  plume,  with  an  increase  in 
AOA/ACO  of  2.3  ±  0.85  in  1.4h  (Yokelson  et  al.,  2009). 

However,  in  the  Williams  Fire,  although  light  scattering  in¬ 
creased  significantly  with  time  since  emission,  AOA/ACO2 
dropped  initially,  followed  by  a  slow  increase  after  ~  1.5-2  h. 
The  net  result  was  a  decrease  in  AOA/ACO2  by  a  factor  of 
0.80  ±  0.21  when  assessing  all  samples  up  to  4  h  since  emis¬ 
sion  (Table  2).  In  addition,  the  assumption  of  a  constant  ratio 
between  light  scattering  and  particle  mass  would  imply  that 
AOA/APM2.5  decreased  from  0.69  ±  0.24  to  0.22  ±  0.04  af¬ 
ter  4  h  since  emission.  This  is  inconsistent  with  the  change 
in  OA  normalized  by  the  sum  of  the  species  measured  by 
the  AMS,  which  decreased  from  a  fraction  of  0.94  to  0.82 
over  this  4h  time  period  (where  the  sum  of  the  measured 
AMS  species  should  account  for  approximately  0.8-0. 9  of 
PM2.5  mass).  Thus  it  seems  unlikely  that  the  relationship 
between  light  scattering  and  PM2.5  was  constant  with  aging. 
We  therefore  cannot  determine  changes  in  PM2.5  over  the  4  h 
following  emission  using  a  constant  conversion  factor. 

Next  we  describe  a  plausible  physical  model  to  rational¬ 
ize  both  the  Yucatan  and  Williams  Fire  observations.  It  is 
possible  that  the  increased  light  scattering  in  the  Yucatan 
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Time  Since  Emission  (h) 

Fig.  11.  ALight  scattering/AC02  (green)  and  AOA/ACO2  (blue) 
mass  ratios  are  shown  vs.  time  since  emission  (h).  ER(PM2.5)  (solid 
green  circle)  and  ER(OA)  (solid  blue  circle)  are  shown  with  er¬ 
ror  bars  depicting  the  standard  deviation  in  the  mean  of  the  ten 
source  measurements.  We  show  all  the  source  samples  to  illus¬ 
trate  the  variability  in  the  source  measurements.  The  log-normal 
(Alight  scattering/AC02)  and  second  order  polyfit  (AOA/ACO2) 
trendlines  shown  best  correlate  with  the  data.  Error  bars  shown  on 
downwind  samples  are  an  estimate  of  the  uncertainty  in  individual 
downwind  samples.  ALight  scattering/ACCL  is  shown  converted 
to  the  mass  ratio  APM2.5/ACO2  under  the  assumption  that  the  ra¬ 
tio  between  bSCat  and  PM2.5  is  constant  with  time  since  emission. 
The  assumption  that  the  conversion  between  fosc at  and  PM2.5  is  con¬ 
stant  with  smoke  age  would  imply  an  unrealistic  AOA/APM2.5  of 
<25  %  at  a  plume  age  of  4h.  The  increase  in  light  scattering  dur¬ 
ing  the  first  1.5  h  since  emission  is  probably  due  to  particle  growth 
via  coagulation  and  then  chiefly  by  condensation  soon  thereafter 
(Sect.  3.3). 


plume  was  due  to  growth  in  the  average  particle  diameter 
towards  the  measurement  wavelength.  The  diameter  growth 
could  have  occurred  by  both  condensation  and  coagulation 
and  the  mass  growth  confirmed  by  the  AMS  in  the  Yucatan 
plume  implies  that  condensation  dominated  over  evapora¬ 
tion  and  coagulation.  The  similar  increase  in  light  scat¬ 
tering  and  the  AMS  mass  may  have  been  coincidental.  In 
contrast,  the  increased  scattering  in  the  Williams  Fire  plume 
could  have  been  caused  by  growth  in  average  particle  diam¬ 
eter  mainly  due  to  coagulation  while,  in  contrast  to  the  Yu¬ 
catan  plume,  evaporation  of  organics  initially  exceeded  con¬ 
densation.  A  simple  polydisperse  coagulation  scenario  can 
explain  the  Williams  Fire  observations.  We  start  by  assum¬ 
ing  that,  at  the  source,  there  are  many  more  particles  in  the 
50-100  nm  diameter  range  (which  do  not  contribute  substan¬ 
tially  to  the  measured  scattering)  than  particles  in  the  100- 
300  nm  diameter  range.  Collisions  between  the  small  and 
large  particles  then  decrease  the  small  particle  number  (with¬ 
out  reducing  scattering),  but  increase  the  average  diameter 
of  the  larger  particles  towards  the  measurement  wavelength 
(530  nm  in  this  case),  thus  increasing  scattering  efficiency 
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according  to  Mie  theory.  Coagulation  may  also  account  for 
most  of  the  increase  in  the  fraction  of  thickly  coated  rBC  par¬ 
ticles  we  observed  in  the  Williams  Fire  plume  (discussed  in 
the  next  section).  At  the  same  time,  ~20-30  %  of  the  OA 
could  have  evaporated  causing  the  small  decrease  in  OA  that 
is  observed.  We  note  that  this  explanation  invokes  the  pres¬ 
ence  of  large  quantities  of  small  (~50-100nm)  particles  in 
the  fresh  smoke.  Particles  ~ 3 0-50  nm  in  diameter  dominated 
the  number  distribution  in  some  laboratory  BB  studies  (Hos- 
seini  et  al.,  2010;  Chakrabarty  et  al.,  2006),  but  there  are  no 
airborne  measurements  of  particles  in  this  size  range  in  BB 
plumes  to  our  knowledge. 

We  also  note  that  the  probability  of  observing  fast,  initial 
coagulation  was  much  higher  in  the  Williams  Fire  plume. 
The  Yucatan  plume  was  from  a  small  fire  and  the  freshest 
sample  was  at  an  altitude  of  1700ma.g.l.  when  the  smoke 
was  already  10-30  min  old  and  the  ACO  reading  on  a  fast 
1-s  instrument  was  823  ppb.  We  did  not  have  1-s  CO  instru¬ 
mentation  during  the  Williams  Fire,  but  our  2-s  CO2  instru¬ 
ment  observed  a  maximum  ACO2  of  390  ppm.  The  ACO2 
maximum  multiplied  by  the  fire-average  ACO/ACO2  ratio 
(0.0713,  Table  2)  for  the  Williams  Fire  implies  a  maximum 
ACO  in  our  freshest  smoke  of  at  least  28  ppm.  Thus,  the 
smoke  at  the  beginning  of  the  aging  sequence  observed  in 
the  Williams  Fire  plume  was  about  34  times  more  concen¬ 
trated  than  at  the  beginning  of  the  observations  for  the  Yu¬ 
catan  plume  and  coagulation  would  likely  have  been  about 
1000  times  faster. 

We  point  out  that  the  decrease  in  AOA/ACO2  in  the 
Williams  Fire  plume  continues  only  to  the  samples  col¬ 
lected  around  ~  1 .5-2  h  following  emission.  After  this  time, 
AOA/ACO2  increases  at  about  the  same  rate  as  scattering 
suggesting  that  some  of  the  organics  that  initially  evaporated 
from  the  particles  may  have  undergone  homogeneous  oxida¬ 
tion  in  the  gas-phase  leading  to  products  with  lower  vapor 
pressure  that  then  recondensed  on  the  particles  (Robinson  et 
al.,  2007).  In  any  case,  condensation  is  evidently  dominating 
over  evaporation  from  1 .5-2  h  onward. 

The  OA  evolution  observed  in  the  Yucatan  plume  is  more 
likely  to  be  the  norm  for  most  global  biomass  burning  than 
the  OA  evolution  observed  in  the  Williams  Fire  plume  since 
the  Yucatan  observations  were  made  in  the  tropical  bound¬ 
ary  layer  where  most  of  the  global  smoke  evolution  occurs. 
However,  the  difference  in  OA  evolution  between  the  Yu¬ 
catan  and  Williams  Fire  plume  presents  an  opportunity  to 
better  understand  the  factors  controlling  OA  evolution  in 
smoke  plumes.  Here,  we  propose  that  the  different  balance 
between  evaporation  and  condensation  in  the  Yucatan  plume 
and  the  Williams  Fire  plume  could  be  rooted  partly  in  the  re¬ 
spective  ambient  conditions.  In  the  case  of  the  Williams  Fire, 
the  smoke  penetrated  above  the  boundary  layer  and  diluted 
in  the  free  troposphere  where  there  was  low  RH  (~1 1-26  %, 
Table  1).  In  the  Yucatan,  the  plume  diluted  in  the  bound¬ 
ary  layer  at  70-100  %RH,  where  particles  would  be  more 
likely  to  have  thick  aqueous  coatings.  Hennigan  et  al.  (2008) 
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found  that  particle  water  soluble  organic  carbon  (WSOC,  a 
surrogate  for  SOA)  may  substantially  increase  upon  reach¬ 
ing  a  threshold  RH  of  ~70%  (though  particle  liquid  water 
content  may  be  strongly  affected  by  particle  hygroscopicity 
and  not  just  RH).  Furthermore,  aqueous  phase  SOA  forma¬ 
tion  results  in  a  higher  OA/OC  ratio  than  straight  gas-phase 
SOA  formation  and  the  high  OA/OC  in  the  Yucatan  plume  is 
consistent  with  some  of  the  SOA  being  formed  in  the  aque¬ 
ous  phase  (Lim  et  al.,  2010;  Yokelson  et  al.,  2009).  The  con¬ 
centration  of  the  hydroxyl  radical  in  the  plume  should  also 
affect  how  fast  the  gas-phase  organic  precursors  for  SOA  are 
oxidized.  In  the  Yucatan  plume,  the  measured  OH  concentra¬ 
tion  was  nearly  twice  our  estimate  of  the  OH  in  the  Williams 
Fire  plume.  Both  higher  OH  and  the  higher  light  intensity  ex¬ 
pected  for  the  tropical  (Yucatan)  plume  should  speed  up  SOA 
formation  in  both  the  gas  phase  and  aqueous  phase.  This 
speculation  is  further  supported  when  comparing  the  rate  of 
ozone  formation:  AO3/ACO  increased  ~4  times  faster  in  the 
Yucatan  plume  than  in  the  Williams  Fire  plume  (Sect.  3.2.1). 
Thus,  conditions  in  the  Yucatan  were  significantly  more  fa¬ 
vorable  for  the  formation  of  SOA. 

We  have  not  ruled  out  the  possibility  that  the  different  OA 
evolution  observed  between  the  Williams  Fire  and  Yucatan 
plumes  could  be  partly  due  to  the  different  fuels  and/or  differ¬ 
ent  volatilities  of  the  initial  POA.  In  a  laboratory  study  of  ag¬ 
ing  BB  smoke  Hennigan  et  al.  (2011)  observed  the  most  SOA 
formation  (growth  factor  of  almost  three)  for  black  spruce 
fuel  and  minimal  SOA  formation  for  white  spruce  fuel.  Thus, 
nominally  similar  fuels  sometimes  gave  widely  different  out¬ 
comes  in  an  experiment  where  the  environmental  variables 
were  monitored,  but  not  fixed  to  a  narrow  range.  In  assess¬ 
ing  global  SOA  production  from  biomass  burning  it  may  be 
significant  that  the  Yucatan  plume  evolved  under  conditions 
with  more  light,  more  OH,  and  more  chemical  pathways  for 
SOA  formation  than  the  Williams  Fire  plume  and  substan¬ 
tially  more  SOA  was  formed.  Finally,  we  point  out  that  the 
large  increase  in  light  scattering  with  plume  age  observed  in 
both  the  Williams  Fire  and  Yucatan  plumes  has  important  im¬ 
plications  for  visibility  impacts  and  climate  forcing  regard¬ 
less  of  the  underlying  mechanism  (McMeeking  et  al.,  2006; 
Reid  et  al.,  2005a).  More  research  would  be  useful  to  deter¬ 
mine  if  the  optical  properties  of  moderately  aged  smoke  are 
better  suited  for  assessing  the  visibility  and  climate  impacts 
of  biomass  burning. 

3.4  Evolution  of  the  black  carbon  mixing  state 

A  fraction  of  0.54  ±0.07  of  the  initially  emitted  rBC  par¬ 
ticles  are  “thickly  coated”  as  shown  in  Fig.  12  (solid  cir¬ 
cle,  see  Table  2).  The  remaining  particles  are  classified  as 
either  “thinly  coated”  or  “not  coated”  at  the  time  of  emis¬ 
sion.  However,  the  fraction  of  thickly  coated  (fTC)  particles 
grew  quickly  and  a  fraction  of  ~0.70±0.10  of  rBC  parti¬ 
cles  were  thickly  coated  within  the  first  hour  following  emis¬ 
sion.  The  fraction  of  thickly  coated  rBC  particles  for  both 
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Fig.  12.  Fraction  of  “thickly  coated”  (fTC)  rBC  particles  vs.  time 
since  emission  (h).  Individual  source  sample  fTC  fractions  are  not 
shown  and  we  instead  show  the  mean  fraction  of  source  samples  at 
t  =  Oh  (solid  circle)  with  an  error  bar  reflecting  the  standard  devi¬ 
ation  of  the  mean.  The  downwind  error  bar  is  an  estimate  of  the 
uncertainty  in  the  oldest  individual  downwind  sample.  Points  were 
fit  to  a  logarithmic  trendline. 


the  fresh  and  1  h  old  emissions  was  higher  than  in  SP2  mea¬ 
surements  from  urban  sources  and  the  rate  of  increase  in  the 
Williams  Fire  plume  was  also  faster  than  was  observed  in 
those  urban  plumes  (McMeeking  et  al.,  2011b).  Schwarz  et 
al.  (2008)  also  found  a  higher  fraction  of  rBC  particles  were 
thickly  coated  in  BB  plumes  when  compared  to  urban  plumes 
that  they  measured  nearby.  Kondo  et  al.  (2011)  used  an  al¬ 
ternative  analysis  to  show  that  rBC  coatings  in  boreal  BB 
plumes  increased  over  several  days,  in  some  cases  leading  to 
increases  in  the  volume  of  the  coatings  by  a  factor  of  two. 

Similar  coating  development  most  likely  also  occurs  on 
the  nascent  inorganic  and  organic  particles.  The  coating  for¬ 
mation  on  the  rBC  particles  (a  flaming  product)  likely  formed 
in  large  part  due  to  coagulation  with  OA  particles  produced 
by  smoldering  combustion.  In  addition,  since  ammonium 
and  nitrate  exhibit  rapid  initial  increases,  the  coatings  likely 
contain  some  inorganic  material  from  gas-particle  conversion 
involving  precursors  produced  both  by  flaming  (NOx)  and 
smoldering  (NH3).  Kondo  et  al.  (2011)  speculated  that  OA 
could  be  redistributed  from  non-rBC-containing  particles  as 
the  plume  aged  through  evaporation  and  recondensation  and 
they  also  showed  that  coagulation  alone  can  lead  to  substan¬ 
tial  increases  in  rBC  coatings. 

Regardless  of  the  exact  mechanism,  as  the  rBC  (and  other) 
particles  become  coated  they  may  scatter  light  more  effi¬ 
ciently  due  mostly  to  size  increases,  but  also  due  to  changes 
in  the  index  of  refraction  of  the  particle  surface  layer.  Thus, 
the  observed  increase  in  rBC  coatings  is  consistent  with 
the  observed  increase  in  Alight  scattering/AC02  (Sect.  3.3). 
The  evolution  of  the  BC  mixing  state  is  also  important 
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because  the  addition  of  organic  coatings  can  increase  the 
mass  absorption  efficiency  of  the  BC  particles,  though  the 
magnitude  of  this  effect  may  be  lower  if  the  coatings  ab¬ 
sorb  light  (Lack  and  Cappa,  2010).  In  addition,  coating  for¬ 
mation  with  salts  and  soluble  organics  increases  the  hygro- 
scopicity  of  the  particles  (Petters  et  al.,  2009;  McMeeking  et 
al.,  2011a;  Zhang  et  al.,  2008)  and  this  likely  increases  their 
contribution  to  indirect  effects  on  clouds  (Koch  et  al.,  2011) 
as  well  as  decreasing  their  lifetime  and  availability  for  long 
range  transport  to  sensitive  snow/ice  covered  regions  (Petters 
et  al.,  2009).  Finally,  black  carbon  affects  human  health  and 
visibility,  and  contributes  both  to  warming  due  to  its  abil¬ 
ity  to  absorb  light  and  decrease  surface  albedo  (Ramanathan 
and  Carmichael,  2008),  and  cooling  due  to  indirect  effects  on 
clouds  (Koch  et  al.,  2011).  Our  observations  in  the  Williams 
Fire  plume  suggest  that  the  evolution  of  the  BC  mixing  state 
could  significantly  affect  the  overall  impact  of  black  carbon 
emitted  by  biomass  fires. 

3.5  Dispersion  of  the  low  altitude  plume 

Our  measurements  of  the  high-altitude  (1200-1900  m 
a.m.s.l.)  plume  generated  by  the  Williams  Fire  detail  the  evo¬ 
lution  of  the  bulk  of  the  emissions  produced  by  the  fire  and 
these  emissions  would  contribute  to  the  burden  on  the  re¬ 
gional  air  shed  along  with  many  other  sources.  As  discussed 
earlier  (Sect.  2.2);  a  small  fraction  of  the  smoke  from  the 
Williams  Fire  detrained  below  600  m a.m.s.l.,  drifted  south¬ 
east,  and  settled  into  low  lying  areas  (Fig.  2).  We  observed 
that  this  smoke  was  noticeably  whiter  in  color  and  was  proba¬ 
bly  more  dominated  by  smoldering  emissions,  which  implies 
a  significant  difference  in  chemical  composition  (Bertschi  et 
al.,  2003).  Flight  safety  considerations  prevented  us  from 
sampling  this  low-lying  smoke.  Thus,  we  examined  data 
collected  by  county  health  departments  (http://www.arb.ca. 
gov/aqmis2/aqdselect.php)  for  possible  signatures  of  the  low 
altitude  Williams  Fire  plume.  The  nearest  air  quality  sta¬ 
tion  southeast  of  the  Williams  Fire  was  in  Santa  Ynez,  but 
only  O3  was  monitored  there  and  no  clear  influence  of  the 
Williams  Fire  was  evident.  The  next  closest  station  was  in 
Santa  Barbara  where  PMjo,  PM2.5,  and  O3  were  monitored. 
A  possible  influence  of  the  Williams  Fire  low-altitude  plume 
was  evident  only  in  the  nighttime  PM2.5  maxima;  which 
were  20,  29,  and  13pgm-3  for  16,  17,  and  18  November, 
respectively.  The  modest  increase  in  PM2.5  the  night  after 
the  Williams  Fire  (17  November)  may  have  been  partly  at¬ 
tributable  to  the  fire,  but  the  peak  was  below  the  air  quality 
standard.  In  contrast,  an  extreme  peak  in  PM10  (1 18  pg  m-3) 
on  the  night  of  18  November  that  was  not  accompanied 
by  an  increase  in  PM2.5  was  almost  certainly  due  to  non¬ 
combustion  sources.  A  comprehensive  network  of  ground- 
based  samplers  would  have  been  needed  to  quantify  the  dis¬ 
persion  and  chemistry  of  the  low-altitude  smoke  and  deter¬ 
mine  the  contribution  of  the  Williams  Fire  to  pollution  lev¬ 
els  on  the  ground.  The  expected  post-emission  transport  of 
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smoke  is  already  a  critical  factor  when  choosing  prescribed 
burn  locations  and  times  in  the  US  (Hardy  et  al.,  2001).  This 
case  study  illustrates  that  air  quality  predictions  could  benefit 
from  research  that  elucidates  both  the  rate  of  smoke  produc¬ 
tion  and  the  distribution  of  injection  altitudes  as  a  function  of 
time. 


4  Conclusions 

In  this  work  we  describe  detailed  pseudo-Lagrangian  mea¬ 
surements  of  the  chemical  evolution  of  an  isolated  biomass 
burning  plume  generated  by  a  prescribed  fire  in  a  temperate 
chaparral  ecosystem.  Despite  occurring  approximately  one 
month  before  the  winter  solstice,  the  plume  was  photochem- 
ically  active  and  significant  amounts  of  ozone  formed  within 
a  few  hours.  Rapid  formation  of  organic  acids  was  also  ob¬ 
served  with  the  precursors  likely  being  unidentified  NMOC 
species.  The  measured  decay  in  ethylene  was  consistent  with 
an  in-plume  average  OH  of  5.27  (±0.97)  x  106  moleccm-3. 
This  is  a  factor  of  two  or  more  lower  OH  than  some  obser¬ 
vations  in  tropical  BB  plumes,  but  sufficiently  high  to  af¬ 
fect  large  changes  in  plume  chemistry  within  several  hours. 
The  calculated  loss  of  NOx  due  to  the  reaction  of  NO2  with 
OH  (~54  %)  and  the  observed  formation  of  PAN  (22  ±  6  % 
of  initial  NOx)  accounted  for  all  the  observed  NOx  loss 
(80  ±  10%)  within  the  measurement  uncertainty.  Roughly 
half  of  the  gas-phase  NH3  was  also  lost  within  the  4  h  since 
emission.  While  both  particle  ammonium  and  nitrate  in¬ 
creased  rapidly  in  the  aging  plume,  the  amount  of  the  in¬ 
crease  was  significantly  smaller  than  that  required  to  balance 
the  loss  of  the  associated  gas-phase  precursors.  We  specu¬ 
late  that  much  of  the  “missing”  ammonium  and  nitrate  could 
have  been  present  in  small  particles  that  grew  slowly  due  to 
the  low  RH  ( 1 1-26  %)  and  were  thus  too  small  to  be  detected 
efficiently  by  the  AMS.  In  addition,  a  charge  balance  sug¬ 
gests  that  perhaps  ~40  %  of  the  ammonium  was  pairing  with 
organic  anions  that  were  also  not  measured.  In  contrast  to  the 
rapid  initial  growth  in  some  particle  inorganic  species,  evap¬ 
oration  of  OA  dominated  condensation  for  the  first  1.5-2  h 
since  emission.  Thereafter  OA  increased  slowly,  but  overall 
AOA/ACO2  decreased  by  a  factor  of  0.80  ±  0.21  over  4  h.  In 
contrast,  in  an  isolated  tropical  BB  plume,  rapid  SOA  forma¬ 
tion  was  observed.  We  speculate  that  drivers  for  the  higher 
SOA  formation  in  the  tropical  plume  include:  higher  OH, 
higher  incident  UV  radiation,  and  higher  RH  (70-100%) 
where  the  latter  would  promote  more  SOA  formation  path¬ 
ways.  The  sum  of  the  excess  values  of  the  species  measured 
by  the  AMS  normalized  to  excess  CO2  decreased  slightly 
over  4h,  suggesting  a  net  loss  of  particulate  matter.  How¬ 
ever,  Alight  scattering/AC02  increased  by  a  factor  of  ~2.50 
over  4  h,  which  may  have  been  due  to  mostly  to  rapid,  initial 
coagulation-induced  increases  in  particle  diameter.  The  in¬ 
crease  in  scattering  has  implications  for  climate  forcing  and 
visibility  even  if  it  accompanied  a  decrease  in  particulate 
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matter.  We  found  that  as  much  as  85  %  of  the  BC  parti¬ 
cles  may  have  been  thickly  coated  after  4  h  following  emis¬ 
sion;  most  likely  with  a  mix  of  organic  and  inorganic  species. 
The  coatings  are  likely  to  strongly  influence  the  BC  radiative 
forcing,  increase  the  BC  hygroscopicity  and  cloud  impacts, 
and  reduce  the  BC  lifetime. 

Supplementary  material  related  to  this 
article  is  available  online  at: 
http://www.atmos-chem-phys.net/12/1397/2012/ 
acp-12-1397-2012-supplement.pdf. 
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